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FOREWORD 


This Manual was prepared by the Structural Research and Development Group, Structures 
Section, Research and Development Division of Republic Aviation Corporation. The work was 
initiated under Contract AF33(616)-6066 in the 750A Applied Research Program, the Mechanics 
of Flight, Project No. 1367, Structural Design Criteria, and Task No. 14002, "Structural 
Analysis Methods". The work is now documented under Task 136710. This work was initiated 
under the direction of the Structural Analysis Unit, Structures Branch, Aircraft Laboratory, 
Directorate of Laboratories, Wright Air Development Center.* Mr. I. Winnegrad acted 
initially as project enginecr and was succeeded by Mr. C. Richard. The Manual was com- 
pleted under the direction of the Structural Analysis Unit, Configuration Research Section, 
Structures Branch, Flight Dynamics Laboratory, Deputy Commander/ Technology, Aero- 
nautical Systems Division, with Mr. G. E, Maddux as Project Engineer. 


The work was coordinated and supervised by Dr. R. S. Levy, Head of the Structural Re- 
search and Development Group, His valvable suggestions and criticisms are gratefully 
acknowledged as are those of the following personnel of the Applied Research and Develop- 
ment Division of Republic Aviation Corporation; Mr. A. Alberi, Acting Manager of Tech- 
nical Engineering; Mr. C. Rosenkranz, Acting Chief Structures Engineer; and Mr. C, 
Meissner, Principal Structures Engineer. 


NOTES ON USING THE MANUAL 


This Manual consists of five basic sections, divided into numbered sub-sections 
and paragraphs, For simplicity in cross-referencing material in the text, all portions of 
the Manual designated with a two-tier number (e.g., 1.1) are considered sub-sections, and 
all portions designated by numbers of three cr more tiers (e.p., 1.1.1 or 1.1.1.4) are 
considered paragraphs, 


Throughout the Manual, the numbered paragraphs (or sub-soections) have been used 
as the basis for numbering figures, tables, and equations, with new sequences beginning 
with each numbered paragraph, Figure and table numbers consist of an appropriate para- 
graph number, followed by a sequence number for the particular figure or table, For 
convenience the paragraph designations have been omitted from the equation numbers, When 
an equation from another paragraph is cited in the text, the number of the paragraph in which 
that equation occurs is also cited. When a paragraph number is not given in conjunction with 
the citation of an equation, it is to be assumed that the equation is included in the paragraph 
in which the citation occurs, 


References are listed at the end of those gections which have more than one reference, 
In addition, eacn section contains its own complete table of contents and list of symbols, 


*Now under direction of Flight Dynamies Laboratory, Directorate of Acromechanics, 
Acronautical Systems Division. 
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ABSTRACT 


This| second volume of tke Thermo-Structural Analysis Manual considers 
additional problems in te field of thermal and mechanical stress analysigjnot 
fully treated in Volume I. Speciallem phasis #6 given to nonlinear analysis. of 
beams and plates and to axisymmetric thermo-elastic analysis of thin shells 4 
Following the format of Volume I, nondimensional graphs, formulas and tables 
are developed where feasible. For clarification of the analytical techniques and 
the use of the numerical data, illustrative examples are given. 


Thd following problems age treated in five individual sections: 


(1) \Large deflection analysis of straight elastic beams with axial end 
restraint and axial end loads coupled with transverse leading and 
temperature, || 

sa ” 

(2) | Approximate determination of tke axial end loads and deformations 

in heated beam columns with initial eccentricities ped 


(3) Y- Approximate solutions for tke buckling of eccentric columns accom- 
“modating nonlinear stress -~strain lawsy , ‘ 


ws 


(4) \ , Axisymmetric large deflections, of circular plates subjected to 
“thermal and mechanic al loads ,.. 


a 


(5) } Axisymmetric thermo-elastic analysis of thin shells. i 


PUBLICATION REVIEW 


This publication has been reviewed and is approved, 


FOR THE COMMANDER: 


ey =: se 


BR. Fe HORNER 
Chief, Strustwas Branch 
Fiight Dynamics Leboratary 
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INTRODUCTION 


eee, 


This Volume TI of the Manual is an extension of the basic work presented in 


Volume Tand covers the following special problem areas in heated beams, plates and 
shetls: 


(4) Nonlinear beam analyses including the effects of initial eccentricities and 
nonlinear material properties. 


(2) Elastic circular plates loaded and heated axisymmetrically, considering 
the effects of large deflections 


(3) ‘The axisymmetric thermoclastic analysis of shells developed in Volume I 
is gencralized, removing the restrictions on geometry and including the 
effects of arbitrary temperature variations through the thickness. 


The material is presented in the form of five independent reports or sections. 


The emphasis has been placed on the development of analytical techniques and 
formulas with their corresponding physical interpretations. Where fcasible, the in- 
vestigations are self-contained, starting with fundamental theoretical considerations 
and notions in the field of static thermo-structural analysis. However, the reader may 


find it useful to refer to Volume I where a systematic development of the basic concepts 
is ‘siven, 


Brief s.mmaries of the five sections of this volume of the Manual follow. 


Secticr 1. Beam Columns Subjected to Elevated Temperature and Mechanical Loads 


This section treats beam columns with axial end restrain - ..:: + nd loads coupled with 
transverse loads and temperature gradients. Numerical results, «+ ot" within the frame- 
werk of large deflection beam theory are tabulated for certain case: 


Section 2. Approximate Solution for an Axially Restrained Column Subjected to 
Elevated Temperature and Lateral Load 


This section presents an approximate methoc of solution of the heated beam colum: 
problem. T!:r method is amenable to problems involving spanwise variations of load the rmal 
gradients an’: stiffness and permits treatment of initial eccentricities. 


Sectior 3, :spproximate Solution for the Buckling of Eccentric Columns 


This s ‘ction develops and presents nondimensional curves to , edict column buckling 
I: rds. Non! near stress-strain relationships are accommodated as well as lateral loads and 
ir ‘tial eccen ricities or thermally induced deformations. 


Section 4. Axisymmetric Large Deflections of Circular Plates Subjected to 
Thermal and Mechanical Loads 


dne m.c..7*9n of membrane stresses and bending is considered for axially restrained 
ce’ rcular plates subjected to heat and load. The differential equations gov 2rning the axisym- 
: atric case are derived and an iterative digital schere is used to obtain numerical results 
‘ra special case over a wide range of temperature and load parameters. 


nt 
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Section 5. Axisymmetric Stresses and Deflections in Shells Due to Thermal and 
Mechanical Loads 


The general equation for linear elastic analysis of axisymmetric shell problems is 
developed including the effects of temperature and load. The special cases of conical and 
cylindrical shells are developed and a numerical example is given for a cyl:. der with 
temperature gradients through the thickness and along the length, 
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BEAM COLUMNS SUBJECTED TO ELEVATED TEMPERATURE 
AND MECHANICAL LOADS 


by 
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Manuscript released by the authors February 1962 for publication as a WADD Technical 
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SECTION 1 
BEAM COLUMNS SUBJECTED TO ELEVATED TEMPERATURE 


AND MECHANICAL LOADS 


1,1, SUMMARY 


This report considers the nonlinear analysis (large deflections) of beams with axial 
end restraints and axial end loads coupled with transverse loading and temperature. 


In order to keep the number of parameters within practical limits the following con- 
ditions are investigated: 


(1) Distributed transverse loads are uniform over the span while concentrated 
loads are at the midsnan. 


(2) Tne temperature varies linearly through the depth and is constant in a span- 
wise direction. 


(3) The beam is assumed to be simply supported at its ends for bending and 
elastically restrained axially (VWigure 1. 3-1), 


Tables cf numerical results in nondimensional form are presented for the cases of 
zero an full axial end restraint in rectangular beams. These tables may be used te detey~ 
mine maximum deflections and bending moments. 


1.2 INTRODUCTION 


In structural analysis beam columns differ from simple beams in that the addition of 
axial end loads and restraints interact with transverse loads in a nonlinear manner, thus in- 
validating the principle of superposition (Reference 1-1), 


When the axial end loads are specified and the beam is unrestrained axially, the solu- 
tion for bending moments and deflections are obtained, in general, by solving a linear, non- 
homogeneous differential equation with constant coefficients subject to appropriate boundary 
conditions, If, in addition, the beam is restrained axially, the total end load is an unknown 
and an additional compatibility relation must be employed, 


The analysis presented considers a beam of constant cross section for which the 
Bernoulli-Euler assumption of classical beam theory is employed. This implies that plane 
sections perpendicular to the centroidal axis before bending remain plane and perpendicular 
to the deflected centroidal axis. it is further assumed that the material behavior is linearly 
elastic. 


It can be shown that the seemingly approximate technique used below (Subsection 1. 3) 
yields results which are identical to those obtained by a moderately large deflectic . analysis 
in which nonlinear strain-displacement relations are employed. 


1.2 (Conttd) 
The following symbols are used throughout this section: 


Width of rectangular beam 

Bepth of beam 

Spanwise coordinate 

Vertical deflection 
Nondimensional central deflection 


Cross sectional area 

Young's modulus 

Magnitude of axial load in beam 
Moment of inertia 

Spring stiffnes.: of axial end restraint 
Half beam length 

Bending moment 

Nondimensional central bending moment 
nown applied axial end load 
Concentrated midspan transverse load 


IZ VUZ SCR aS “i x oo 


Q,Q Nondimensional concentrated midspan load 
T,, TS Temperatures at lower and upper beam faces, respectively 
Ty Ty Nondimensional temperature differences between upper anc lower beam faces 
a Nondimensional average temperature 
WwW Intensity of uniformly distributed transverse load 
W,W Nondimensional uniformly distr.buted load 
rag Coetficient cf linear thermal expansion 
B Ratio of distance from lower beam face to centroidal axis to the total depth 
i 

+ an 
a isl 
X \L, nondimensional 
g ot Tess 

SUBSCRIPTS 

M Due to mechanical loads 
T Due to temperature 


1,3 DERIVATION OF BASIC EQUATIONS 


The bear (Figure 1.3- 1) is referred to a rectangular coordinate system in which 
transverse deflections y and distances along the centroidal axis are meusured from a point on 
the undeflected centroidal axis at mid-span. Specified axial end loads are denoted by P and 
the elastic axial restraints are shown schematically as springs with stiffuess 2K. We assume 
uniform distributed loads of intensity Wand a mid-span concentrated load of magnitude 2Q. 
The temperature variation through the thickness ts linear, varying tram YT, aethe lower ex- 
treme fiber to a at the upper extreme fiber. 

Since the remperature distribution varies linearly with respect to Cartesian voordi- 
nates it produces no stresses in an externally unrestrained beam (Reterence 1-2). ta this 


cage, the stress-free thermal curvature is piven by 


fee Oe Q) 
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FIGURE 1.3-1 BEAM COLUMN - MODEL USED FOR ANALYSIS 


We now consider that the external specified and redundant mechanical loads are ap-~ 
plied to this thermally bent beam, Designating the additional curvatures preduved by these 
mechanical loads as Ye it follows that 


ee \ one & Ld ao iu 
YM “EL” BI E Hy - (WL + Q) (i =x) F (L- J , (2) 


where positive moments cause compression in {he outer niin y is the total deflection, and 
H is the magnitude of thy axial load in beam*, The negative and positive signs in the first 
term of the right hand side of Eq. (2) refer to compresgion and tension respectively, Adding 
Eqs. (1) and (2) and rearranging yields 


Hy . Ww 2 (WL EQ (=x) Op | 
OR Co ampere an El nite 1p 4 (2) 


Due to symmetry it {sg only necessary to consider one half of the beam, so that the dilferen- 
tial equation given by Eq. (3) applies tn the Interval Qs x ~ Land y(x) is an even function, 
i,e., ¥(x) = y(-x). The boundary conditions are given by 


y' (y= yiby= 0, (4) 


1+ SOLUTION OF EQUATIONS 


ee ne at eee nent ne eretarr me ene rt ern 


The complete lutions of Eqs. (3) and (4) of Subsection 1.3 for the cases of axial 
compression and tension, respectively, are written as follows, 


* His always taken as positive, regardless of whether the axial load is tensile or compres 
sive 
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(1) Axial Compression 
ase Soa eg oe 
EI) dcosdA L EIR 
cia) 
2 r ‘ 2 
secre a “3 a5 Ese Se eet eee ie 
2EIX EIX A : EIA 
(2) Axial Tension 
ue 2 3 sr a SBE se . (Ty ~ 40) 
EIX A coshAL | ELA 
(ib) 
we Qe leer gy Mh we we 
2EIY" EI? 4 h‘"o i EI 281 ED? 
where 
Az 3 for vompression, 
EI ay 


~,\/H eal evans 
dR: +\ BT for tension, 


The solution is rat yet complete since, in general the axfal end load H and henve A is unknown, 
An additional compatibility relationship must be ~mploved to evaluate this quantity, Such a 
relationship tg obtaiped oy noting that rhe chanye to the spanwise distance between the beam 
ends due to banding, thermal expansion and mechanical axial strains must be equal to the total 
change tn length of the axial end restraints, This condition can be expressed by 


] ; j Gs 
x (v')" dx + nl - Ver 4 (1 - R) 1) eS ee ce (3) 


j 


where positive and aeyative sipas associated with H refer to comipression and tension respec- 
tively and 8 defines the location cf the centroidal axis (Figure l.4-1), 


centroidal axis 7 


eee 


PIGURE T4-1 PARAMETER fo, DEFINING CENTROIDAL AXIS LOCATION 


be 


Ld 


1.4 {Conttd) 
Substitution of Eqs. (1a) and (1b) into (3) and. simplifying yields for the compression case 
& freind) (cos? 4 - 8 + 3] 
U oy a 
r 
x 2 
‘Es (ais 3 
y = a =, 
2 4X Beets x3 x2 X : 
~~ ee | vs 
Se a) ee 
: x L A cosy ss 
Ow in} x 1 
+ ae : eons cos. : 3) (4a) 
A r IN 
( —_ - |W . Qsinx Ty 
= - (sink ~ A cos A) [mt =) rapa 
cos i 
a 
, ie : mH Po Ps oe io ee 1 | 
AL? KL 
and for the vase of axial tension 
2 j = 9 Jy S ay | 
ace in X) (cosh™ A. 2) + | 
\ a 
{AL sinb2at } ot | | W_ Qsinka , a? 
2 AX ae ay. 3 Xx 
eosh™ 7X 4 
$2 ) 0 
oe, | Beton. 2] (.%  Qsimn | 
<4 | <3 <3 <2 \ 
3X ~ Yo cushy X ) (4b) 
oan eer /sinb 3 cosh Ya \ 
og 4 = + 7 2 “oof 
r | A eee & 
ow rs ee We Q ainh \ t. 1 
pS eI ACCOR A Se eS eee 
\eosi: 4 ye os 
ew " 
° ~h : ah) Vos ee 2 [#1 +(l- 8) T | 
2 a | s { 
AL Ki. KL. 
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where 


2 
ee SG 
AL= + vj for tension 


=e —— for compression 
EI 
aL (T, : T,) 7 % 
h ae oe 


Equations (4a) and (4b) are transcendental equations from which \ may be determined, 
Equatiens (la) and (1b) then yield the deflections for the known values of 4. In general, kqs. 
(4) are difficult tc solve except by graphical or trial and error procedures, In addition, pre- 
sentation of numerical results for the large number of parameters appearing in this veneral 
problem is cumbersome, For these reasons, the number of parameters te be used in the pre- 
sentation of numerical results have been reduced by restricting considerations to the following: 


(1) The beam cross section is rectangular and therefore pe > 

(2) The beam ts prevented from moving axially at its ends (K -* ~ ) or free to 
Move axially (K= 0), 

(3) The transverse load is either uniform over the span or concentrated at the 


mid-spat., 


1.5 NUMERICAL RESULTS FOR BEAMS OF RECTANGULAR CROSS SECTION 


CASE A - Uniform transverse load over the span with ends rigidly restrained axially: 


The beam is shown schematically in Figure 1.5-L. For this case, we substitute 


rar 
\ 


| ] : , nae ; . , ; 
=u, Rt, and K =O into Eqs. (lL) and «iy of Suleseetion 1.4.01 cing rondimensional 


quantities: 


a iS 


ae ty Noa 


Roa 
~ PHL : : sak 
\ 4 Vi ae (tension positive) 


1.5 (Cont'd) 


FIGURE 1,5-1 RECTANGULAR BEAM WITH UNIFORM LOAD OVER 
THE SPAN AND ENDS RIGIDLY RESTRAINED AXIALLY 


(1) 


ss 
tt 
| rca | 
=< 
; a | 
i 


| d-e2k | | 1 | [3-2 i. we 
2 4X 4 Sena) te d J ai 
(2a) 

a = isin \ - X cog X) EB . 7,| + ae =T 

vos A ce 
=. A Ww W 
ne E = | Ee | <2 
: (3a) 
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for the case of compressive end loads, while for the tensile case 


1 sinh 2) 1 Wy, wee 1 Ww 
sis br aioe i 7a me ae, ee 
i J LA“ cosh*A d 3X 
(2b) 
= = =2 
=e = q = as 
+ —?¥— (sinh X- dooan) |Z + T,| -AL=7 
X° cosh % N 
uae eh Wy, 1 W 
‘hecaet e jeer eae Soc) ay 
52 | 32 d coshyv 9X2 
[ (3b) 
<7 Wow 1 
M = - |-¥ +5 ee ait) 
{ x2 d | cosh A 
Val of X, y, and M for various combinations of Ty T and W are given in Table 1,5-1, 


This table permits the determination, either directly or by interpolation, of the axial end 
loads, maximum (central) deflection and maximum bending moment corresponding to speci- 
fied temperatures and transverse igading, A typical case, extracted from the table is plot- 
ted in Figure 1,5-2, showing the variations in end load, central bending moment and deflec- 
tion with average temperature when the temperature difference and transverse load are 
held constant, 


The figure shows that extremely large values of T are required to raise the compres- 
sive end load to values in the neihborhood of the critical Euler value, AoR zs 
\ ~ 
This is due to the fact that additional beam expansions caused by increasing the average tem- 
perature are accommodated by further bending with very little change of compressive end 
load, 


ae Ont 


Thus, although the theoretical end load in a perfectly steaight axially restrained 


: , es a ; 
column reaches the buler buckling load when T = 4 ofl, values of T several orders 


of magnitude higher than this can actually be achieved without excessive stresses or de- 
flections vecuring, tn Chis respect thermal buckling differs significantly from buehlinge 
caused by deadweight fonds, since stresses and deflections tend to become excessively large 
when mochanieal loads exceed the Euler buckling vahie. 
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FIGURE 1.5-2 NONDIMENSIONAL END LOAD, CENTRAL DEFLECTION 
AND MOMENT VS, NONDIMENSIONAL AVERAGE TEMPERATURES 
Wear, F212 
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CASE B - Concentrated transverse load at midspan with ends rigidly restrained axially 


For this case, we substitute W= 0, & a and = =0 into Eqs, (1) and (4) of Sub-section 
1.4. Nondimensionai quantities are as defined by Eq. (1) and in addition we define 
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for the case of compressive end toads, while for the tensile case 
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TABLE 1. 5-2 
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(Pages 1.38 through 1,55) 
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VALUES OF X, y AND M FOR SPECIFIED VALUES OF T 
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CASE C - Transverse load uniformly distributed or concentrated at midspan - speci~ 
fied axial loads (zero axial end restraint): 


The beam is shown schematically in Figure 1,5-3, Since the ends are unrestrained 
axially, the beam ends are free to move due to the action of temperature, transverse loads 
and specified end leads. Hence the compatibility Eqs. (2) and (5) do not apply, and since \ 
is now a known quantity, the maxirnum deflections and bending moments may be decermined 
directly from Eqs. (3) and (6), These quantities are independent “‘ the average temperature 
in the beam. 


2Q |b 


FIGURE 1.5-3. RECTANGULAR BEAM WITH ZERO AXIAL END 
RESTRAINT AND SPECIFIED AXTAL LOADS 


It is important to note that when the end loadings are spacified then a modified princi- 
ple of superposition as shown by Figure 1.5-4a may be used, The figure shows that tho re- 
sultant effect of several transverse loads and temperature distributions acting slmultaneously 


Q ye Q W 
H HoH debe cs ea | H 
b Bh t: y te # t he 


(Temp, distribution Ts") (Temp, distribution "T ") (Temp. distribution 
2 mS + "To") 


(a) Superposition valid for an axially unresgirained beam with specified end loads "H" 
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(Tamp. diatribut{on Ty") (Temp. distribution "To (Temp. distribution 
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(%) Superposition net valid for an axially restrained beam 


FIGURE 1.5-4 MODIFIED SUPERPOSITION PRINCIPLE 
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in the presence of a specified axial load can be obtained by superposing the effects of the indi- 
vidual transverse loads and temperatures acting with the axial load, 


However, this superposition principle is not valid when the ends are restrained axially 
(e.g., Cases A and B), In such cases the axial loads depend on the transverse loads and tem- 
peratures nonlinearly, 


L6 UsE OF THE TABLES 


Tables 1.5~-1 and 1.5-2 are reproductions of IBM 7090 digital computer numerical 
solutions for the beam-column problem, ‘Table 1.5-1 presents results for the case of uniform 
transverse loading and Table 1,5-2 applies for concentrated midspan loads, Each table is 
first subdivided into sections corresponding to given temperature differences and transverse 
loads and then further subdivided into compressive and tensile end loading cases. Since the 
quantities W, Ty and T specified in a given problem will not in general coincide with those 


listed in the tables, interpolation must be employed. The important quantity to be evaluated is 
the nondimensionai end loading parameter, Xand the spacing of the tabulated values has been 
made sufficiently close so as to allow reasonably accurate interpolation for engineering pur- 
poses, A systematic interpolation formula will be given as part vt an illustrative problem 

in Sub-section 1.7, 


Numerical values are given in terms of a floating decimal number system and are to 
be interpreted as shown by the following examples 


0, 2456582E 00 = 0, 2456582 x 10” = 0, 2456582 
0, 2456582E 02 = 0.2456582 x 10” = 24. 56582 


0, 2456562E-01 = 0. 2456582 x io”? = 0, 02456582 


1.7 NUMERICAL EXAMPLES 
EXAMPLE I - Boam with Full Axial End Restraint: 


Figure 1, 7-1 shows a simply supported strip with immovable ends subjected to a uni- 
Ib, . , : 
formly distributed load of Zin . The temperature varies Unearly through the thickness from 
100° F at the upper face to 150° F at the lower face, Young's modulus ar: the Unear coeffici- 
ent of thermal expansion are taken to be 


B= 30x 10° pal 


a =6x io in/in - °F, 


Find the axial end load, midspan deflection, bending mement and the maximum stress. 
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‘ 10" 


100°F 
EZZZTTEZD | ~s. 


aT nes a 150°F 


FIGURE 1,7-1 ILLUSTRATIVE PROBLEM; STRIP WITH IMMOVABLE ENDS 


SOLUTION 


From Figure 1.7-1 and the given data: 


b = yn T =100°F 
oO 
h = 0,10" Ty = 150°F 
ae E =30 x 10° psi 
Yb 
vV- 9 —_ -§ 
Wese2ee in a =6x107 in/in- °F 
Therefore 
4 
— sy. =) =x SEBO: a) ° 
tro = ~ ae «SS 


(30) cos® (1) 
» 


2 2 
Pang L + y -p 4,75 ‘é ~15 ee 5 
Frsa () er-t,) =6 any a) (100-150) = - 0.75 


iH 


=y 2 6 / Ss 2 
Tao (=) cri +t,) =6 oy? (2. 1004150) = 3.75 
. Vi ) (T i ) eee) 0, 10 noe ee 
Using the above nondimensfonal quantities we must now utilze Table 1,5-1 to determine the 
nondimensfonal axial loading parameter A. Since the table does not vive e Addit ectly for the 
above combination of quantities the following interpolation procedure {gs recommended: 


(1) Determine the next Jowey and higher va'ues of W that are given in the table, 
Designate these values as Wo and Wy respectively. For this example 
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(2) Determine the next lower and higher values of Ty that are given in the table, 
Designate these values as T ao 20d TT ay Tespectively, 
For thic Usample 
Ts =-4,0 and 
Da Oi 
(3) For each of the four combinations (Wo, Tad), (Wo tas) (Wy, a), 


(Wy. Top the table lists X's corresponding to Tts on either side of the given 


value for Tf, Denote the values of T lied in the tables for the combination 


(W,, Ty) by (Ty, Th) and the corresponding values of \ by Ou Mi) where 
TT tt. 
Y i 
For this example 
‘F  =3 39 Xt _ 
Co0 3.397 An V2 
ry = Yt a 
rho 3, 896 rO0 =1,1 
Toy = 3,041 Mgr eke 
PT aa oe 
HS a 017 1,3 
Po =3.416 Any = 1,5 
7 me gee yt 
Teas 1O= 14 
Ty, 73-948 Rye 8 
"Tt sd. fd Xt 
TM , =444t Meas 3 
(os) The value of 4 my new he obtained from the following interpolation formula, 

W- © Ww, - WW 

x = ace Soni ae A + quik S367 B 
wo. WW : | Wy - Woy { 

1 0 ‘ee J 
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where 


Substituting the known quantities into the above formula yields 
X =-.524 


Thus the axial end load is given by 


ie 6 3 nT Re Z 
H= oA = (30 (10) C1) (. 10) /125 (~2 524)" = 27,4 1b. (compression) 


L (5)7 


The nondimenstonal deflection and bending moment at midspan can be determined from 
Table .5-1 using an identical interpolation procedure. However, a simpler and more ac- 


curate method is to substitute the known values of X, W, and Ty inte Kas. (3a) of Subsection 
LS*). Direct substitution yields 


rae L = 1,60 and 
i Xe 0 
fem? | 
M = oman 20504 %, 
isbh Pear 


39 that 


= (.10) (1.60) = .16 in. (downward) 


s 4 —_ 
M] = ath M 
eee 12 


‘ 6 4 
= ee (2, 94) 
12 (5) 


= 29.4 in.-lb {compression in upper fiber) 


* Eqs. (3b) are not to be used here sinc. chey apply only for > 6 (tension) 
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The maximum stress is thus 


oa wk + M® ~ 17 900 psi. 


The effect of axially restraining the beam is evidenced by comparing these results 
with those for an axially unrestrained and simply supported beam subjected to the same 
transverse loading and temperature but with H= 0, In such a case the central deflection, 
bending moment and maximum stress are given by (Section 4 of Reference 1-2): 


: 2 
Spwiigt > 2 UMS tg) te 
| _,  S84El 2h 
= (104 + ,038 
= . 142 in 
wet)? 
M] =* = 25 in-lb, 
8 
= O 
(6) = 15,000 psi. 
max 


Thus, in this case, neglecting axial end restraint which may be present yields unconserva- 
tive results for both the central deflection and maximum stress. 


EXAMPLE UW - Unrestrained Beam Cclumn With Prescribed Axial Load: 
Figure 1, ¢-2 shows a simply supperted strip with movable ends and a preseribed 
20 Tb tensile load subjected to a concentrated midspan load of 10 lbs, The femperature varies 


nearly through the thickness from 200° F at the upper face to 15° F at the bottom face, 
Young's modulus and the linear coetfictent of thermal expansion are given as 


(a) ; 

BK= 30x 100 psi 
-O. oo, , 
a=6x 10 in’in- °b 


Find the midspan deflection, bending moment and the maximum stress, 


a8 


Qué8 20? 


Looe 
ua 


BPIGURE L.7-2 TLLUSTRATIVE PROBLEM, STRIP WITH MOVABLE ENDS 
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SOLUTION 


From Figure 1, 7-2 and the given data: 


h= jn T> =200°F 

h = . 08 T, =150°F 

L= 6" 

E =30 x 10° psi 

Q=5 Ib 

H = 20 lb ese e infin Ser 
Therefore 

= 12QL° mes 

Q= 8h. - 2 6) 6 _ 10.55 


Ebh 30 (10) (1) (.08)4 


2 


t 


¢6f6\ 
(T, - T,)=6(20) von i) (200 - 150) = 1,69 


a ee 
a7 8 () 


Es - | 20.19" we 
(30)(10)° (1 ¢, 08)°/12 


ot 


As discussed in Sub-section 1.5, since the ends are free to move, the bending response 
of the beam is independent of the average temperature. The central deflection® and bending 
moment can thus be obtained by direct substitution of the nondimensfonal parameters into 
Eqs. (6b) of Sub-section 1.5. This yielus 


y 3 | =. 674 
zs x 3.0 
Tan 
M lM Le = > OA 
iebh Fic 
so that 
v | =, 054 in (downward) 
Jx=0 
M =26,4 in-Ib (compression in upper fiber), 
-X=0 


1.8 ADDITIONAL CONSIDERATIONS 


(1) Application of Basic Equations to the Case of Unequal Elastic End Restraints. 


The general formulas presented in Sub-section 1.3 have been developed for the case 
of equal elastic end restraints of stiffness 2K (Figure 1.3~1), As is frequently the case, these 
restraints may be unequal, This situation can be accomodated by replacing the quantity 2K 

4K_K 


in Sub-section 1. 3 with the equivalent quantity ana where 2K, and 2K, are the unequal spring 
1 2 
stiffnesses. Thus for example in the special case where one end is held (B, =<) then 
lim 4K,%, 
pe a eS 


= 4K; ana this quantity is to be substituted for 2K in Eqs. (3) and (4) of 


Sub-section 1.4. This procedure will yield the correct results for the deflection and 
bending moment. 


(2) Maximum Bending Moments and Deflections, 


This report presents results for the determination of the central deflection 
and bending moment for a symmetrically loaded heam-column, It can be shown by the usual 
maximum-minimum procedure that the central bending moment is always numerically the 
largest bending moment inthe beam, If the temperatures and loads tend to produce curva- 
tures in the same direction then the central deflection ig a maximum, However where 
temperature and loads tend to relieve each other, the central deflection may not be a maxi- 
mum, i this quantity:is desired it may be found by considering the full deflection formula 
(Eq. {1) of Sub-sectior 1. +4). 
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SECTION 2 


APPROXIMATE SOLUTION FOR «.N AXIALLY RESTRAIN“D COLUMN 


SUBJECTED TO ELEVATED TEMPERATURES AND LATERAL LOADS 


2.14 SUMMARY 


Approximate solutions are obtained for the axiai compressive load and the lateral 
deflections of columns which are restrained bv an axial spring and subjected to an increase 
in temperature, The solutions are presented graphically in terms of nondimensional para- 
meters. 


The preblem considered is a column of arbitrary cross section with pinned or clamped 
ends subjected to an arbitrary temperature and lateral load distribution in addition to flaite 
initial eccentricities. The effects of the thermal pradients, lateral loads, non-linear axial 
springs, and plasticity of the material are discussed, The simpier problem of constant 
bending stiffness is explored to illustrate the evaluation of the nondimensional parameters, 


2.1.1 Definition of Symbols 


The following symbols are used throughout this section: 


¢ Amplitude of lateral deflection, inches 
b Axial extension parameter (b = aT ¢ 
b Width of rectanzular cross section, inches 
i A, \ 
d Axial shortening parameter (42) | 
\ i fp: 
h Depid of cross section of eclumn, inches 
i Integer indicating order of the deformation mode 
K Axial stiffness cf end restraint, pounds per inch 
t Leneth of column, inches 
in. Coefficients expressing Che taveral deformation as a polynomial, inehes 
j : } 
q Lateral lead, pounds per inch 
t foree funerien, Katie of caial load in column toa reference 


load (r(x) = Pex} FO) 


Ww Lateral deflection of column, inches 

x Axial coordinate ef colunim, inches 

7 Lateral coordinate af coltumep cross section, mehes 

i Coefficient of thermai Gineur) expansion, inches per ineh. °F 
ay Bellection of emits of cetunin, inches 
A 
€ 


Incremenial change 
Axial strain due io axial load or tersperature, inches per ined 


: : : An ‘ z 

€; Buckling strain correspendimus ty i mode; X ro) \ 
* Load parameter ( n ‘ \ 

ee ee 

\ 
” Curyature of column, 1b inehes 
A; Ligenvalue for which non trivial solutions of tne differential equilibrium 
: BN 
aqlation exist i 1 } , 1/ésq. inches 
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Figenvalue (H; ~ tan H.) 


a 


Nondimensional axial coordinate (€ = x/£) 


Effective radius of gyration = \EVEA , inches 
Axial stress (¢ =F/A), psi 
oO Stiffness function, Ratio of bending stiffness of cross section te a reference 


vending stiffness (< (x) == ts) } 
‘oo 


QU Mm EF 


S 


é \ 2 
Additional axiai shortening function (i = Gar oot 


, ‘ .d, 

oe ast ki ere Te PO Gee 

¢ Nondimensional linear axial shortening term | @=1+ == +2 — 
\ x aor) 


A Cross sectional area, sq. inches 

C Amplitude of curvature (v= re, %) 

ES Secant modulus ae are i fora linear material), psi 

pee, ‘ ae i < U \ “ . 

EA Axial stiffmess (EA = 2 E dA, note (Eo zdA=), Ib 
$s S 

ae . oe race! y Ms 

BI Bending stiffness (EP = 2 Bo 27 dN), pounds square inches 
eos 

b Axtal toad in column, tb : 


cary = . « a “t 
Moment of inertia of cross sections, Inebos 
Hait length of column (Lh. = 2), inches 


\ Moment (4b+ Bbw, pounds per inch 
p Redundant transverse ead, ib 
iG Temperature inerement from datum, *F 


SUBSCRIPTS 


re) Condition before appl cation ef axial load 
i a mode 

d Due to lateral mechanical bouds 

- Die to comperatiure 

Sus, Siraple simple bousdary conditions 

Soe, Simypoe clumped pe Cary cenditions 
Clu. Clamped clumped Geo sidiary conditions 

s Svermetrical mede 

a AnGiesy nimepricab mode 

Bo GENERAL APPROACH 


fotoriratien: character stics 


The analysis et the problem is depron hed Dy oxaniing Une : . 
of the column with a linear material to oobtain . sefaton which sat. feos ceubilitn and 


equulibriam, Phe tellowin. steps ure emp over to selye the colum show Fisare “bod, 


oe) The column detormes in characteristic iades weich are depemtent upor the das 
tribution of the * endime stitfmess and axial load inthe cubona as owell as the boundary gend) con 
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ditions. Basic orthogonality relationships exist among th characteristic deformation modes 
which are used In obtaining a solution for the column in a thermal and mechanical environment, 


(2) The deformaticns of a coiumn subjected to lateral loadings and temperature 
gradients in addition to an axial load are then determined by solving the equilibrium equation. 
Yhe anulysis results in a solution in terms of the "initial'' deformations caused by the lateral 
loads, thermal gradients, and manufacturing and loading eccentricities acting on a column 
with no axial load. These "initial" deformations grow with the magnitude of the axial load. 
Each characteristic deformation mode is magnified as a rectangular hyperbola, by a factor 


( “ which is a function of the applied load and the characteristic "buckling" load cor- 
NN 7 * / 

responding to the deformation mede, The method of determining deformation modes and cor- 
responding buckling loads is illustrated for columns of constant bending stiffness with pinned 


or clamped ends. 


(3) The axial load is the: xpressed as the solution of a compatibility equation 
which considers the thermal expansion, the lateral deflection, and the axial deformation of « 
restraining sping, Tne axial lead determined by this method satisfies both the equilibrium 
and compatibility conditions of the structure. The solution is, therefure, correct for a re- 
versible (one to one relationship of load and deformation) structure even though a column de- 
forms non-linearly with load. 


(4) Tne compatibility equation is quite difficult to solve but a simple and reason- 
able approximation of the additional axial shortening of the column due to the axial load re- 
duces the compatibility equation to a form which car be readily solved by graphical means. 


(5) The solutions of the axial load and deformation of columns are presented for 
pinned (simple) or clamped ends simplified by applicable formulae and graphs together with 
a computational procedure and illustrative problems. The techniques to obtain the "initial" 
deformations due to lateral load and thermal gradients, snd the effective stiffness of the 
restraining spring are also presented, The effects of nonlinearity in the spring or material 
are discussed. 


2.2 ANALYSIS 


The mathematical relationships necessary to solve the problem are derived in 
Reference 2-1, The results are summarized below, 


2,2.1 Relationships of Deformation Modes 


Various orthogonality relationships exist between the characteristic deflection modes 
and their derivatives whenever the boundary conditions are natural ‘>, g,, free, clamped, 
simple, etc.). These relationships are useful in evaluating the initial and final lateral de- 
formations as well as the axial shortening in termis of these characteristic modes. 


Orthogonality 1 glationships, (see Section 1A of Reference 2-1) are obtained by solving 
the homogeneous differ: tial equilibrium equation 


(own) + Arwl)t = 0 (1) 
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where 
rv) is the ratio of bending stiffness to a reference bending stiffness (EQf)) 
WwW is the lateral deflection 
r is the ratio of the axial load in the cross section to the reference 
load (Fy 
Fo 
r is the eigenvalue EY . 
eo 


The following orthogonality relationships are derived 


4 
| (rwi)! w, dx = 0 ixfk (2) 
oO 


Thus the characteristic modes, wich are the solutions to the differential equilibrium 


equation, are orthogonal to the derivative of the weighted slope foe wi)! | . For end loads 

(r = 1), the deflection and curvature modes are orthogonal, This relationship permits the de- 
termination of the amplitudes of the deflection modes for the lateral deflections caused by the 
lateral ioad and thermal gradients and permits the determination of the growth of these modes 
when the column is compressed, 


L 
I, rw! whe dx =0 ik (3) 


Similarly, the slopes of the deflection modes are orthogonal witn respect to a load 
weighting factor (r). This relationship permits the rapid determination of the axial shortening 
of the column in terms of the magnitude of the deformation modes. 


a 
owt wt dx = 0 ixk (4) 
fe) ik 


The orthogonality of the cu -vature with respect to the stiffness weighting function «@), 
permits the solution of the non-homogeneous differential equilibrium equation resulting trom 
the initial deformations due to the laterul load and thermal gradients, The effect of the 
luteral load and temperature can be expressed as initial curvatures which grow as a rectang- 
ular hyperbola with an increase in axial load as indicated in Figure 2, 2,1-1, 


The homogeneous differential equilibrium equation can be rewritten as a function 
of the curvature (x =w") 


ile. (ex trArx=0 (1) 


The solutions of this equaiion results in eigenvectors of the survature which are or 
thogonal with respect to the weighting function @ (Eg. (4) of Parag: aph 2.2.1). it is assumed 
that these curvature modes form a closed set for the "natural" boundary conditions (columns 
of constant ET @p =1) result ip Fourfer expansions). Thus, any curvature can be exp: essed 
as a weighted sum of the eigenvectors. 
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i.e. c= eh X (2a) 
where (ox,)" + h; rx, = 0 (2b) 


If the unrestrained column is not straight because of initial eccentricities, thermal 
gradients, and la..ral loads, then the moment acting on the column is proportional to the 
change in curvature and the equilibrium equation hz:omos, for an end loaded column (r = 1), 


Cr HE iv ( a Xe) - (3) 


where 2) is the curyaluve caused by the iatera! load and Map is the curya- 
ture vareed py ihe thermal gradient, 

Expressing the solution as a weighted sum of the characteristic curvatures (Eq. (2a)) 
which satisfy the homogeneous equation and manipulating Eq. (3) results in the solution for the 
curvature of an axially loaded column in terms of the initial curvature and the load ratio. 


Bi 
X= I=, x; (4) 


wth . 
where be corresponds to the i eomponent of the initial curvature of the column due to the 


lateral load and temperature. Thus the column defornis under axial load by having each com- 
ponent of the initial curvature increase by a magnification factor 1/(1-7,). Correspondingly 


the slopes and laters! deflection modes increase by this same factor, 


Since the initial deformations are expressible as an infinite series of the character- 
istic curvatures, it becomes expeditious to employ the initial curvatures and to integrate 
this infinlte series to obtain the expression for the slope which is employed in obtaining the 
axial shortening. Employing the initial slope or lateral deflection will result in a siower 
convergence of the series expressing the axi.l shortening and would reduce tue accuracy of 
the approximate procedure employed to solve the compatibility equation. 


2,2.3 Determination of Deformation Modes and Buckling Loads 


The deflection modes (w) of the column are obtained vy solving the homopencous 
equilibrium Eq. (1) of Parapraph 2.2.1. The buckling loads (Fi) are those values of the axial 
load which result in non-trivial solutions of the equation, The technique of obtaining these 


iteras is illustrated for a column of corstant bending stiffness ( =1). 


The homogeneous equilibrium equation for constant bending stiffness (E}) ts 
w 
Ses EN Beye BS oye) (3) 


The general solution is 


WES og Ry P Op cos f/X xX + C4 sin Ax (2) 


| ie, 
where Cyr Cyr Cgy Cg are constants to bs determined by the boundary conditions, 
2,8 
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Solutions for various boundary conditions can be found in various texts. 


simple, clamped-clamped, and simple-clamped are summarized below. In addition the 
(1) 


The simvle- 
method of determining the magnitude of the initial curvature modes is indicated. 


The boundary conditions are 
w(0) = w(t) =0 


w'"(0) = w"(z) =0 


‘ sacar : — Pe 
“.W,= a, sin =a, sin Pu ee (3a) 
her = : 
where w yy, 
KF. . 2 
on ieee (it ‘ 
El r, \ t ) (4a) 
. tmx ss 
4 = sin ares (oa) 
er 
x. dx = f/2 
J 1 
0 
ad 
| wit dx ; 
; “0 ae . imKx 
|e ee ew Sin Se x (6a) 
j i: L 
t » tg 4 
| uF dx 
0 
where w" 2. Cm 


(2) 


Clamped- Clamped Column 


VYhe boundary conditions are 
w(0) = w(d) = 0 


wi(Q) = wh) = 0 
This results in two different types of modes. Modes which are syminetrical about the 
mid-length of the column and modes wh'ch are anti-symmetrical. 


For symmetrical modes 


Meee 


2iuax 
a, (1 - cos -——-—~) 
is i 


(3b) 
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eae a: ° n (4b) 
El “i 2 
= 2imx 
Xe cos z (5b) 
L 
[ wit %G ax ‘ 
iC: =~ _=5/ w'' COs RELL (613) 
is 9 Lido L 
* dx 
0 


For anti-symmetrical modes 


vu X- sin an. & 
au - sin au 


x 
1L- cos 2u. ~— 
a aoe it i 
ii 2p, ~ sin 2u, 1 - cos 2u, 
x x 
sin 21“ cas 2h L 
Me poo Sse = tan 
i 2, ~sin Zu, i cos Zu; a 


An alternate form, employing the mid-length of the co amas an orivin, results ina 
simpler result which is identical to the joining of wo simple-clamped columns, 


~ (3e) 
as oe 
sin > 


» ‘ 
: ja oi ia, ao ae (4c) 
kil i y 12 

where tan Be =U; 


ant th 
and y= 1.43 u,- (2 > ) 
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Hi 
Mea sin o a (5c) 
| 4/2 wx, dx rl ‘ 
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paw 2 Oakes 
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(3) Simple Clamped Column 


The boundary conditions are 
w(0) = w"(0)= 0 


w(t) = wi) = 0 


/sin s 
ie oe (3a) 
im“ \sing, 2 ‘ 


where tan Ue = 


2 

OF u; 

ho 22 (4d) 
. x = 

x sin Wim: (Sd) 


3 o{ . 
Cot ae Bed saath aces a wwe ) ow" sin u, park 
i 2 98 2. y+ - 
Vee ace. aie Mil ery oe (6d) 
“0 2 ! 


The deformation nodes and buckling toads are required inthe determination of the 
end load and deflection of a restraiaed column as indicated in Subsection 2.3, lo the boundary 
conditions or stiffness destribution do not conform to the case summarized above, then the 
buckling loads and modes must be determined for the coluann cunsidered, Special cases may 
be found in various text and recourse must be taken fo approximate sohitions when the column 
becomes more involved, 


The method of solution for a column with a linear materiaf, propo: -doin Chis report, is 
quite general and can handle the more complex contigurations provided the initial deflection can 


be dese ribed by characteristic eigenvectors with known cigenvalues, 


2.2.4 Compatibility Equation 


The axiatly unrestrained column (Figure 2. 1-1b) will dotorm when subjected to lateral 
loads and temperature. This wiil cause an axial motion of the eads of the column of Bes 
ayy 


. 
bo =tat- A, (1a) 
where aT = € T is the unrestrained axial strain due to temperature. Ht the strains are not 
uniform then 
1 pt Jorvda 
a T= z i) coe, —— dX (ib) 
“0 | dA 
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A = axial shortening of the column due to the lateral load and temperature 
© gradients through the cross sections 


~ 


at f° we! +wit)“d (1c) 
o 2 JO T q x 


The application of the axial load causes the column to shorten due to the axial strain 


and additional axial shortening (A - A). The final movement of the ends must be equal 
to the deformation of the axial spring F/k, 
© AG Sh es ee ee (2a) 
7% o 8 6AE a) k 
Dividing by the quantity te) results in 
aT Ao Pr A 4, FE 
Res ie nan Cee ee ee 2 
ey AE ey €; k €} : 2b) 
i sy ae ae — X 2 %O) 
Wate Sy AR A = Oe (2c) 
Noting that 
F oa bea € : 
reas = Es == Are) SS ST) (3a) 
AE €) AEE, Ak & Cy 1 
and 
Fo Uo KeA _ ELA . 
Rte, kley eT hy) 
and letting 
aT-A f& 
ie) i 
ee b (4) 
1 
and evaluating the axial shortening 
= 1 4 f t 1 2 = 4 4 t 2 as 1 if t 2 ta 
A.=5 is Guia oat de J. we dx = 5 ' (wt) dx (Sa) 


but because of the orthogonality of the characteristic slopes (fq. (3) of Paragraph 2, 2.1) we 
obtain 


a 


d.=5 . w 2 axed s f wr sta 5b) 
o 2 eae saa ge des io? O% 9 
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and 
: 2 
Be ee 1 pt 8 1 | 
= = ee t x t = f -—_—— = 
A-B,= 5 [> [twn?- wh y?] ax=> [oor — L| ax 
fo) a J 
(5c) 
1 ae 
Sous eee ra 1 
= > fe ) 1 i (wi) 1x 
letting 
1 rd 9 
= | ere 
A oh do MI) ey 
ee ed (6a) 
and 
pe are 
17, ames (6b) 
we obtain 
4, 
——— ve 
te wwe) 
and 
A- 4, 
te, nd Si (6d) 
Substituting Eqs. (a), (8b), (4), and (6 into Bq. (2b) results in 
KA = “ 
b - ny - rdie, id ie pag (7) 


Tho above equation is quite difficult to solve, but the solution can be simplified if it is 
noted that the expression de, denotes the lncrease in axial deformation due to the g -owth of 


the i mode of deformation. The first niode is magnified to a much greater depree than the 
higher modes when the axia' load is compression, UW fact the first mode becomes predominant 
as the axial load approaches the first couckling load.  Vhe approximate method utilizes this 
fact in approximating the value of ado, hy expandiys dio, (for ic: 1) as a power series and 


only employing the most significant term and by a. proximating the infinite series by a finite 
series, This approximation mei ‘dis not recommended for tension loads, since the approxi- 
mation for oO and the infinite series may be incerrect, 
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Noting that 
4 1 2 2 
el CO i ie, (8a) 
and that 1 ae eas {i = 2) 
© © n (8b) 
; ey qs @; =d,9, * Dz die, is approximately dO, + 2 a 7,4; 
= {=2 i=2 
(8c) 
where nis a sufficiently large integer. 
Equation (7) can thus be approximated by 
n 
/ ¥ Yr 
SA | me _ i 
b ~ d,@, - on tee 7, 4) my (9a) 
Letting 
n 3 
EA : ‘i ¢ 
ee inh Tees OR a 9 
¥ + Aaa) my Si (9b) 
{=2 
we obtain 
3 
b Be 
eeoN Ss to Gan es = Ye 
x Cr ae" (te) 
or 
b - day - 1h = 0 (Id) 
where 
i _ z (9e, 
and 
we oe 
d= xk (9 f) 


The value of the axtal strain parameter (n,) can be determined from a graphical plot 


(Figure 2.2,4-1 and -2) of the variation of b with ny for various di This value of n, can then 
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be employed to obtain the axial load acting on ine column and the deformations of the column, 


it should be ncted that the approxirmate method results in an upper bound on the axial load since 
it ignores a small part of the change of axial motion due to the higher modes. 


2.3 PROCEDURE 


The information and technique needed to solve the problem are enumerated in the fol- 
lowing sections. 


2.3.1 Determination of Nondimensional Parameters 


By the definition Eq, (6a) of Paragraph 2, 2. 4: 


Equation (9,1, 2 - 3c) of Reference 2-2 indicates that 
row? a feo? oe 
Bh. ee ee a ak 


aes ogee a apse: 
a ae 20, oy 


: wt)? dx 


14,3 + (1b) 


and 


(2a) 


Noting that 


1 
— Sa SS pe eee (3) 
id 
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m9; 
d and ar will be evaluated for the chosen boundary conditions and for constant El (@=1) 
° 1 
utilizing Eqs. (5) and (6) of Paragraph 2, 2,3, 
(1) Simple -Simple Column 
5 
o? 1/2 a accra 
i — ¢ i 
cS oe ae tee (2b) 
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9 4 n 
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re th as aan ra 
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(3) Clamped-Simple Colurmn 


ats O08 u; 4 Je acs 
i 2 bi. i 2 iL 
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(4) Summary 


Equation (9a) of Paragraph 2.2,4 c#n then be summarized for columns of con- 
stant EY as follows: 
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at -A,/4 Ces EA 4 eo ‘ 
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2.3.2 Computational Procedure 


The following computational procedure is recommended to obtain the approximate 
solution for the axial load and lateral deflection of a restrained column. 


(1) Determine the axial expansion of the column due to temperature and the 
lateral load, assuming, at this step, that the restraint against actual expansicn is removed. 


, 4 
taT-d =] ° eqdx - | axe 
° i £ e wi) dx (1) 
where wo= Ww tw, 
te) q 7 
Nas lateral deflection due to lateral loads (usually obtained from 


reference texts) 


T= lateral detlection due to thermal yradients (method of determin- 
ing Wap indicated in Section + of Reference 2-2), 


€, ae ‘ ee 
T= axial expansion due to temperature (method of determining ¢€ 


: 
indicated in Section 4 of Reference 2-4), 
(2) Calculate the critical axial strain from Eq. (2c) of Paragraph 2. 2. 4) 
\ i 2 
a>. SLAs TLE 
ed 
rw 2. : : 
(@.£. a ay & for simple-simple. go 1) 
(3) Expand the {nitlal curvature (wh ; re, x), due to the lateral load and tem- 


perature, in terms of the characteristic curvatures of the column. A sufficient number of 
terms of the series should be (aken to ensure accuracy. The coefficients of the character- 
istic curvatures are obtained by the general equation 


r¢ 
J} Sy dx 
‘ oC 
a ane (2) 
| Oy dx 
BS) 


Formulas for constant EI (¢ 1) are presented In Eqs. (6a), (6b), (6¢), and (6d) of Paragraph 
2,2.3. Tables to evaluate the inteyrals of Eqs. (6a) and (6b) of Paragraph 2, 2.3 are pre- 
sented in Tables 2.3.3.1, 2-1 and -2, 


(4) Evaluate the pertinent parameters 
.27 2 
te C; jOx, dx 
(a) qd, ig Op (Reterance Eq. (2a) of Paragraph 2.3. 1) 
1 2p AW 
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Appropriate formulas for qd, are presented for columns of constant EI and for the 


chogen boundary conditions in Paragraph 2. 3,1. 


ae 
(a.g. d,= -+- for simple-simple column), 
i 42 : 
4m p i 
n 
EA »: ni 
(b) @= 1+ cet 2 . qd, (Reference Eq, (9b) of Paragraph 2.2.4) 
i=2 ‘1 


Values of € are indicated in Eq. (5) of Paragraph 2.3.1 for constant EL 


4 4 C2 
a EA i x. i ‘ : 
(eg. €=1 + er ace laa for simple-simple column) 
oT 0 =) i 
(c) qd ae (Reference Eq. (9f) of Paragraph 2.4%. 4) 


atT-A at-4 ,, 
= oft oft : : i 
(d) b= mee -“— (Reference Eq, (96) of Paragraph 2, 2,4) 


2 
€,@ 4,0 ® 
(5) Enter Figure 2.2,4-1 or -2 with the appropriate values of b and d, and de- 
termine 1: 
(6) The axial load ig then determined from Eq, (3a) of Paragraph 2, 2.4, 
€ e 2 
F=o0AskE — €,A- EAA, 0D 4 (3) 
€) l 1 1 
(7) The lateral deflection of the column could likewise be determined as 
- vw. 
vt ‘ 
w= : L n, Gate) 
r 
where the value of U Pas ot vgs (413) 
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ee £4 : ‘ oS ips 
and the v>lue of in om be calculated or determined with the aid of Figure 2.2. 1-1, 


i 
reo : 
a ae {ax/! 
(e.g. Ww ee) ee ') AUP ae i for pinued end column). 
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(8) A slightly better approximation to the value of n, can be obtained by employing 


aT - Aa/t 4, 2 
2 é a ° a qd; 
r, p £49 i=l 
ss gy ee (5) 


which corrects for the error in approximating the initial shortening by a finite series. The 
additional shortening of the higher modes is still approximated by the first term of the series 
in order to simplify the solution of the compatibility equation. 


(9) The graphs can also be employed to determine the approximate average temp- 
erature rise (with no cross-sectional gradient) required to produce a maximum permissible 


r 
ste oh, the . . ee : max 
lateral deflection, The permissible ratio of maximum to initial lateral deflection Ce 
) 


cai be emploved with Ficure 2,2, 1-1 to estimate a value of " which can be utilized with dy 


in Figures 2. 2,4 to determine a value of b which could then be employed to calculate aT, (if 
the higher modes of the initial deflection are significant, then they must be included as in- 
dicated in Eq. (4a)), The initial eccentricity will permit the calculation of ® and qd). As- 
suming a value of ny or of aT will permit the determination of the other from Figure 2.2.4. 
Tho value of un could then be used to determine w which could be plotted against aT; aT 


could also be plotted against the axial load by employing Bq. (8a) of Peragre * 2,.2,4. A 
typical plot is shown in Figure 2, 3. 2-1, 


F = F(a) - 


FIGURE 2,3.2-1 AXIAL LOAD AND LATERAL DEFLECTION 
VS AVERAGE TEMPERATURE RISE 
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2.3.3 Initial Deformations 


The initial deformations before the end load is applied must be determined to a fair 
degree of accuracy, Two methods are presented herein. The first procedure assumes that 
an analytical expression for the initial lateral deflections is available as a polynomial, while 
the second method assumes that the lateral deflection is known only at a finite number of 
points (e,g., determined experimentally or by computations), 


2.3.3.1. Polvnomial Representation 


The lateral deformation of a beam can be determined by integrating the equilibrium 
equation (Eq. (3) of Paragraph 2, 2.2) with A =0 and employing the given boundary conditions. 
Solutions for lateral loading are available in literature for the type of boundary conditions con- 


4 3 4 
: - 44 (X)_ 9 f% x a ith unie 
sidered (e.g., Wa 548i G ) 2 ( ) + G ) for a simple-simple beam with uni 


form lateral load), The slopes and curvatures can be obtained by differentiating these poly- 
nomial expressions, The lateral deflection due to thermal gradients is not as readily 
available but can be derived quite simply by utilizing the equations of Section 4 of Reference 
z-2, The technique is illustrated for the simple-simple and clamped-clamped beams, 


dog Be Da hd Lateral Deformation of Beams by Temperature Gradients 


Let the thermal gradient through the depth of the beam be expressed as a 
polynomial in the spanwise dimension, 


ont 


aT -atT, j : 
8.8. 0 i AaT 1 x _i j 
: = 2 > sy (F) = a Em, g (1) 


h A h 


where qT; and y, are the lemperatures on the positive and negative side of the column with a 
linear gradient (see Figure 2, 1-1). 


If the structure {gs statically determinate, as ina simple-simple beam, then 
tho thermal gradient represents the total curvature (wt) since no redundant forces are pro- 
duced by the temperature, The slope and deflection cab then be obtained by integrating this 
cuevacure and emploving the given boundary conditions, 
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If the structure is statically indeterminate, as in the clamped-clamped beam, 
then the curvature is affected by the curvature caused by the redundant loads. Equations 
(4.2,2.5 - 7 and -8) of Reference 2-2 can be employed to determine the deformation of the 
clamped beams where subjacted to a thermal gradient. The thermal gradient aT = 

b 


Bu yi? induces redundant loads Py and M, and the following total deformations result: 
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The polynomial definities of curvature can be converted io a elgenvector ea- 
pansion by means of Eqs. (2a) of Paragrapis °. 2.2 and (2) of Paragraph 2,3, 2. 
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For the simple-simple column with constant EI, this results in the simple 
Fourfer gine expansion of the initial curvature. 
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The clamped-clamped column with constant El and symmetrical loading re- 
sults ina Fourfer cosine axpansion. 
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The curvature coefficients for the clamped-simple and the clamped-clamped 
columns with anti-symmetriecal load must be similarly evaluated, 
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Tables 2.3.3,1.2-1 and -2 present the values of the integrals for the simple- 
simple and the symmetrical clamped-clamped columns. 


TABLE 2.3,3.1.2-1 FOURIER EXPANSION OF MONOMIAL FOR SIMPLE-SIMPLE COLUMN 


SG, i) 
0 6366 “0 | alae | 0 12732 0 
1 3183 -. 1591 .1061 | -,07958 06366 | -.05303 
2 1893 ~. 1591 1013 | -. 07958 06262 | -. 05303 
3 1248 -, 1350 .09894 | -.07655 | .06211 -, 05215 
4 ossi4 | -. 1108 scrat | -.07353 06061 | -, 05125 
es 1 | -.oov7e | .osaea [ -.ososs | .o5862_ |. 05011 
oe _-.o7497 | 07489 | -.06560 | .o56z9 | -, 04872 
ens oe 7, 06208 |, 06647 | -. 06099 | 05368 | -.04712 
8 ~ 05280 05889 | -. 05631 05088 | -, 04537 
“9 | .oa683 | -.04503 | 05223 | -.05176 | .os799 |, 04350 
19 | -.o3877 «| .ogeas | -.0a74n | co4sto | -o4i5s 


9  e Ps 
A | SS eee * ees vy , 
Mus Sate of ae dee eS mS, 1) (2a) 
i“ I aaa a) 
(jy 
oh 4 
su, )= | €! sinimede 
pi (Cy aka " 
key cae a ee 
: oe fie H . aca | : ote : o =f 
Sq. i) = Ae seo os ay! ) Seek (2h) 
Zim) es iy (i 7) VJ ny! 


2.3.3,1.2 (Cont'd) 


TABLE 2, 3.3.1.2-2 FOURIER EXPANSION OF MONOMIAL FOR CLAMPED-CLAAMPED 
COLUMN i 
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2.3.3.2 Discrete Lateral Deformation 


In some instances the lateral deflection is not known as a continuous function, but 
is determined analytically (see Paragraph 4.2.1 of Reference 2-2) or experimentally at a 
discrete number of points. The deformation of 'n" discrete points can be employed to obtain 
an approximating expansion of 'n' characteristic deflection shapes. The amplitude of the ''n" 
characteristic deflections are determined by matching the displacements at the known points, 


The technique is illustrated for a simple-simple column with known deflections at the 1/6 
points (n= 5). 
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The symmetry of the odd characteristic deflections and the anti-svmmetry of the 
even characteristics can be employed to reduce the order of the simultaneous equations 
necessary to solve for the amplitudes (a,). 
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it should be noted that C, = - a, (2) 


i 


Similar solutions can be obtained for different boundary conditions and for dif- 
ferent types of known deformations. 


A rapid solution can be obtained if one assumes that the lateral deflection is pri- 
marily in the fundamental mode (e.g., uniform thermal gradient on a simple-simple column). 
Thus the values of the nondimensional parameters can be simply determined with a minimum 
amount of computation. As an example, the simple-simple column equation can be expressed 
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Similar expressions can be derived for the other boundary conditions, 


2.3.3.3 Axial Shortening 

The initial axial shortening due to the application of temperature and lateral load 
is determined as a function of the polynomial coefficients expressing the lateral deflection, 
The axial shortening fs determined as a function of the actual lateral deflection rather than 
by a function of the amplitude of the characteristic modes, in order tv improve the accuracy 


y 


of the solution (see Eq. 5 of Paragraph 2.3.2). 


As:ame that the Lateral deflection (Wo) is readily available and is expressible as a 
polynomial of the fourth degree or less, Higher degrees of polynomials can be treated tne 
manner similar ’9 that indicated below. Differentiating the deflection results in the slope 
(Wo') which can oO. squared and Integrated over the length of the column to obtain twice the 
axint shortenlag 
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Typical examples for simple-simple columns are presented. 


For Constant lnermal Gradient G27) 
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It should be noted that any manufacturing or loading eccentricities will augment 
the lateral deflection parameters (Cj and dj) by causing an axial shortening when the column 
is loaded, but should not be included directly in the Qo (initial axial shortening because of 
thermal gradients and lateral loads) of the compatibility equation, The ‘sternl eccentrici- 
ties due to the various causes a: idditional and are included in the additional axial motion 
of the ends. The axial shortening however, is not linear but is proportional to the square of 
the deformations, Thus the initial axial shortening (Ao) must be computed as the difference 
of the axial shortening between the column including the manufactured eccentricities and 
the column containing only the manufactured eccentricities (Wao). 


2.3.4 Spring Constant (k) 


The value of k to be employed in the solution of the problem can significantly 
effect the magnitude of the resulting axial load. A column with zero axial restraint (k - 0) 
doe not develop any axial load. The load increases with increasing axial restraint until 
it reaches a maximum for complete restraint (k = ©), 


The value of k is the value of the load produced by moving one end of the column 
a unit distance in the axial direction ruiative to the other end (assuming that the column 
offers no resistance), This is the stiffness coefficient defined in Section 4.2.5 of Reference 
2-2, Tf both ends of the column are spring mounted, then the vatue of k is determined by 


putting the two end springs in series 
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The method of determining the value of k {na composite structure {ts (lustrated 
below for e truss joint. 
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E, = A,E, = AE) 


Yor symmetrical geometry 8, = 6, = 6; Ay 
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For unsymmetrical geometry 


The load will cause a non-vertical motion of center strut and change the angles. 
For small deformations the following relatiunships result: 
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2.3.5 kitfect of Non~ Linearity in Spring or Material 


Solui.cns to the axial load in the column can be obtained even when the axial re- 
straint and/or stiffness of the column varies with the axial load. 

The solution of a column with a variable axial restraint ean be obtained vy super- 
imposing a plot of the flexibility parameter [sare k = EA/Uk (my) | as o function of the axial 
load parameter (71) upon a plot of the solution of the compatibility equation for the load para- 
meter [ny os niCGe) | as a function of the flexibility parameter. ‘The first plot is obtained 


directly from the spring characteristics whereas the second pict is obtained by varying the 
yalue of EA/kt to obtain different values of 41. The intersection of the two plots, as illus - 
trated in Figure 2.3.5-1, will result in a compatible solution. 
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FIGURE 2,3,5-1 FLEXIBILITY VS AXIAL LOAD 
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FIGURE 2,3,5-2 AXIAL LOA VS AVERAGE THERMAL EXPANSION 
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The solution of a column with a load dependent material stiffness (Eg) due to 
the plasticity of the material can similarly be obtained by a trial and error procedure. 
Assuming a value of 71 will determine a value of Eg. This will determine a value of the 

EA 

flexibility parameter ee which varies linearly with Eg; values of the lateral deflection 
parameters (C;) whose lateral load component varies inversely with E, and whose thermal 
component is unaffected; and the axial shortening parameters (dj) which varies as the square 
of the Cj. Corresponding to these values of 7) , values of dj can be computed and em- 
ployed in Figures 2.2.4 to obtain values of b which can be utilized to calculate aT. The 
actua! value of 7)1 resulis when the calculated value of aT corresponds to the correct 


value of aT, A graphical approach is illustrated in Figure 2.3. 5-2, 


This above method cr be utilized to solve a problem with a load~dependent 
axial restraint and a load dependent material stiffness. The above technique assumes that 
the axial and bending stiffness of the column is not significantly affected by the variation of 
the stress through the thickness of the cross section but is determined primarily by the 
mean stress (I/A), 


2.3.6 Illustrative Problem 


The computation procedure for the approximate solution is relatively simple 
and is illustrated for a problem which can be directly compared to an "exact" solution ob- 
tained from Table 1.5-1. 
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AXIAL EXPANSION 
From Eq. (ib) of Paragraph 2.2.4 


Tht t me 
at =a 25+ = 6(10)° 323 = . 0021780 


AXIAL SHORTENING 


From Eq. (2d) of Paragraph 2.3.3.3 
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EXPANSION OF CURVATURE 


From Eq. (2c) of Paragraph 2.3.3, 3 
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Direct sub: titution in the Fourier sine expansion or the use of Table 2.3.3.1. 2-1 
will enable the determinotion of the magnitude of the curvatures. Direct substitution is em- 
ployed in tnis jllustration to obtain general equations for the case of uniform load and thermal 
gradiencs, Tue Fourier coefficients defined by Eq. (2b) of Paragraph 2.3.3,1,2 are 
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From kg. (2a) of Paragraph 2,3,3,.1.2 
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The values of C; are then calculated. The same values would be obtained from Table 
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PERTINENT PARAMETERS 


¥rom Eqs, (2b) and (4b) of Paragraph 2,3, 1 
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en a, = d. + 4, +d, + dy = $4,513 
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Employing Bq. (9a) of Paragraph 2,24 and including correction for finite sum approximation 
of the initial shortening (Eq. (5) of Paragraph 2.3.2) results in 
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Substituting the calculated values results in 


0 = (264.8 ~ 94.555) - (94,555 ~ 94.513) ~ 92, 895 0, 7 (1. 234) ny 


ao 
i 


= 170,203 ~ 92, 895¢, - (1,234) n, 


137,93 ~ 75.239, - n, = b- dye, - ny 


cn) 
i 


Entering Figure 2.2.4-2 with b = 137.93 and d, = 75,28 r ‘ults in the solution My = 407, 
AXIAL LOAD 


From Eq. (3) of Paragraph 2.3.2 
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From Eq. (4a) of Paragraph 2. 3,2 
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The geometry and loads employed in this illustrative problem results in para- 
meters which correspond to an exact solution found in Table 1.5-1, 


These parameters are 
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<I 
] 


4 
W/L 
= Eb =) = 3.0 


D 
= L 
a a(=) (T,*+T,) = 108.9 


A comparison of the exact and approximate solutions is presented below. 


Exact Solution Approximate Solution 
N - 1,000 ~ 1,005 
14 405 407 
F 4,014 4,024 
Wise 92" oT 


Note that the approximate solution results in slightly larger values, The appro..i- 
mate method, however can be applied to more general types of loadings and boundary condi- 
tions. 


It should be notad that the illustrated example shown in Figure 1, 7-1 of Saction 1 is 
solved by an interpolation of the "exact" tabular values. The approximate graphical pro- 
cedure as well as the interpolation of tables procedure is subjected to relatively greater 
errors in the vicinity of low load ratios, This is because of the greater relative signift- 
cance of the higher modes in the graphic approach and the large slope of ths T vs} {as 
illustrated in Figure 1. 5-2 of Sect? . 1 of the tabular solution which permits large vart- 
ations {nv A for small variations in T. Fortunately the design of the column is not de- 
termined by this condition of low axial load. The interpolation and graphical solutions are 
not subject to significant errors, however, when the ivad ratios are higher. 


A comparison ot the interpolation solution and of the approxlmate graphical 
solution of this report is presented as follows; 
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Interpolated Sol. Approx. Graphical Sol, 
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3.1 


SECTION 3 - APPROXIMATE SOLUTION FOR THE BUCKLING OF 
ECCENTRIC COLUMNS 


3.1 SUMMARY 


a E 2/3 
An approximate solution {F, = oS ee eae J5{ 22 1- =t0) is obtained 
1 hs So'o r Eo, 


for the buckling of an eccentric column which includes the effects of the stress level, the 
stress~strain relationship, and the initial eccentricity. 


The analysis technique is presented in a nondimensional form to permit the analyst 
to consider various cross sections, boundary conditions, materials and degrees of eccen- 
tricity. The type of cross section is reflected in Ure I/A ratio which defines the square of 
the radius of gyration (r). The boundary conditions are -‘ontained in the value of Ky? 


which consider the end fixity and length of the column. The material characteristics are re- 
flected in material parameters (Eq; oo 8) which represent the stress~strain relationship. 


The nondimensional analysis graphs are presented for several values of the eccentricity 
“10 
ratio ( 10 which may oecur because of initial waviness, eccentric loading, lateral loads, 


or thermal gradients, Analysis for intermediate vaiues of the eccentricity can be conducted 
by interpolation. 


Illustrative problems are presented to indicate the computation procedure and the effect 
of the initial eccentricity upon the stability of the structure. 


3.1.1 Definition of Symbols 


The following symbols are used throughout this section: 


Cviginal amplitude (eccentricity) of fundamental mode 


10 
Depth of cross section of the column 

L Length of column 
Ba Coefficient expressing initial deformations of the column as a power series 
q Lateral load acting on column 
v Radius of gyration of cross section 
w Lateral deflection of column 
x Axtal coordinate of column 
z' Nistance from reference axis 
z' Distance from reference axis to bending axis _ ¢ . 
ze Distance from reference axis to centraidal axis, oh = jetda! | dA) 
z Distance from bending axis 
A Area of cross section ct j 
CG, 1) Fourier expansion coeffictent (C(j, 1) = ie EE cos? TEdE) 
= aN Initial slope of the ovs ey curve F 

> Secant modulus of the o vs curve (E, = af 

s ; FBS gS So 
Fa, Secant modulus at the bending nrxis 

, : 3a/ 
ae Tangent modulus, slope of the g vs Ci curve {E. = oe ) 
Eno Tangent modulus at the bending axis 
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a 


KA Axtal stiffness of cross section (EA = | Eg d A) 
J 
x . — ff 2 
El Bending stiffness of cross section (EI = J Eo z dA) 
won = as ek] _s= K 
EI Buckling stiffness of cross section | EI = EI +x Soc Fy 
22 
¥ Axial load on structure (compression positive) 
Ro Buckling load of column with zero eccentricity or a linear material 
a E 2h 
= } oe K 1 : 
F Buckling load of column (F, = —S Er 1 -{/6 40 (, oe 
1 1 2 “Sore r/R 
4 So 
. : _K 
F Euler Buckling load (F.. = —-5 EI 
E E 2 oO, 
4, . 
Te Second moment of area {inertia) of cross section ki = i (z' - Zz) as) 
K Stability constant depending upon boundary conditions and the bending stiffness 
/ BP ok 
My) otw NY dx 
: R L 
8 Meaty ee : 2 &x ) 
distribution along the columus (KC Ra ; pear is ieee 
Oe 7 ee Gee, 
vo Lt 
M Moment acting on crogs sections (Positive M causes compression in positive 
fibers), 
: : , af 3 \ 
Q, Third moment of area of cross section (2, 24 (Zt = z') dA} 
: ee ee aes Deere 
Sj, 1) Fourter expansfon coefficient (SG, o Fesingw Ed € 1 
T Temperature rise above datun 
a Coefficient of linear thermal expansion 
a Parameter expressing variation of secant modulus in the cross seetion 
’ pe Ee 
Ne CNY 
woof flew -f 
ee, ae OS 
.) yo 7 
B Nondimensional parameter employed in mathematical definition of stvess- 
strain relationship (Seetion 3 of Reference 8-0) 
¥ Nondimensional parameter expressing the average variation of the tanjcent 
f Mw DY 
modulus in the eross section | + Ee ale taka 
\ vol Va 
Ao VY Difference in thermal expansion 
9) Operator deneting a small variation 
€ Amtal strain of eross seetion 
Axial strain caused by stress (€ € eb} 
£g : Kg ; 
ey Average axial strain eerresponding te an axial lend Y 
f Axial strain at the bending axis 
S Average axial strain corresponding te an axial lead ae 
” Shift of bending axis (n z = hy 
x Curvature 
£ Nondimengional axial coordinate (€ x» 4 
Cc Axial stress of cross section (compression positive) 
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om Axial stress at bending axis 
To Reference stress in nondimensione! stress-strain relationship (Section 3 
of Reference 3-1) 
oO Nondimensional parameter expressing the variation of the bending stiffness 
along the length of the column (G(x) = EL/E | i) 
E 
@ Nondimensional parameter (»- (F ‘A/o : ‘KI A fod a 
° 
ar Fe 226 E\| 2/3 
rei) Nondimensional parameter ©® -(¢ AS 0 e 
= J 
a 
ge) 
SUBSCRIPTS 
M Caused by mechanical load 
T Caused by temperature 
1 Pertaining to the first (fundamental) mode 
oO Pertaining to initial or datum 


3. 1,2 Discussion of the Problem 

The stability of a structure fs a very complex problem, Exact solutions exist only for 
very few special cases. Approximations must be attempted when the olutier is complicated 
by vatiatlons (reductions) of the stiffness with the applied load because of the resulting non- 
linearity of the equation, The reductfon in the stittness is caused by a non-linear structural 
material, The non-linear stress strain reiationship reduc s the stiffness of the structure by 
lowering the secant modulus and shifting the neutral axts, 


Temperature, time, and the eccentricities of the siructure tend te reduee the allow- 
able magnitude of Che applied Tead on the structure,  Elev.ted temperatures reduce the stiff- 
ness of the structure, facrease the non-linearity of the material because of plasticity vt the 

tresses caused by the applted and thermal stresses, and usually increase the cocentricity 
at the structure. The axial forces acting through the eccentricity of the structure imposes 
Moments upon the structure which increase at a zreater rate than the applied load. These 
moments cause a vartation of stresses through the cross sections which can reduce the bend- 
ine stiffness. The larger the eccentricity the earlier the initiation of the non-linearits of the 
material and the lower the stability of the struciure. The bending stiffness can also decrease 
with time because of creep of the material, 


The approach, emploved tn this section, to ea. mine the stability of s column (Figure 
3.1. 2-1) is to determine when the {nerease of external moment acting on the column cross 
sections tend to become greater than the inerease of internai moment that can be generated 
by the cross sections, The axtal load acting on the column at this time is the buckling: load 
and ts a function of the baundary conditions, the bending stiffness of the cross sections, and 
the rate that the bending stiffiess ta decreasing. The bending suffness is a function of the 
vross=section, fhe magnitude of the stresses and thetr distribution and the stress-strain re- 
tationshio., The bending stiffness usually decreases with an increage in the applied load, 
temperature, time, or eccentricity of the column. Thus the column becomes unstable when 
the applied ioad becomes equal to the buckling load which is uetermined by the bending stifi- 
ness distribution and its rate of decrease, This can occur by increasing the applied toad 
and/or decreasing the magnitude of the buckling load. Vhe buckling lead is reduced by the 
vecentrictty, the temperature, and me as indicated previously. 
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View A-A 


FIGURE 3.1,.2-1 ECCENTRIC COLUMN SUBJECTED TO AXIAL LOALSS 


The evaluation of the buckling load requires the deter: 1ination of the stiffness of the 
column and its rate of change. This is complicated by the shifting of the bending axis towards 
the convex side of the cross section because the higher stress and strain icvels on the con- 
cave side reduce the secant modulus, The stiffness relationships, derived in Appendix A of 
Reference 3-2, are employed with the stability criterion described previously to obtain an 
approximate soluticn for the buckling load of the column in terms of the stress and lateral 
deflections, The solution is then approximated in terms of average stress level and initial 
eccentricity which can be readily computed. The inilial eccentricities may result from the 
unloaded shape, eccentric loading, lateral ioads, or thermal gradients. Methods o: com~ 
puting the initial eccentricity due to load and temperature are presented in Section 2 and 
summarized in this section, 


An analogous approach can be employed to transform the loads or stresses to 
strains. Thus the column becomes unstable when the average axial strain becomes equal 
to a buckling strain. This ap_coach is most convenient for a nondimensional presentation 
of the analysis, 


The solution is predicated upon the following assumptions: 
(1) The bending strains are small compared to the axial strains in the column. 


(2) The temperature distribucion can result in deformations and reductions in 
the moauli but do not cavse sipnificant thermal stresses. 


(3) There is a one-to-one relationship between the stress (a) and the strain fe ) 
for each fiber of the column, This aesumes that there is ne stress reversal Chroughou 
the loading history. h is also assumed that the stress-strain relationship can be approxi: 
peas 
mated to a sufficient degree of accuracy by the relationship : as - {1-f) "Fo, 1 Bsinh e oe 
aaa) : 


defined and deseritbed in Section 3 of Reference 3-1, 


The approximate solution, which results from these assumptions, is of a 
form which does not violate known solutions of special cases and is indicative of the struc- 
tural behavior in experiments. 


3.2 ANALYSIS 


A column subjected to loads and temperature (Figure 3. 1.2-1) must satisfy equilib- 
rium, Thus the change in applied moment (Mad must be equal to the change in the internal 


moment (M, The changes in external moment arise from “hanges in the axial load (F) 


int’ 
and the lateral deflection (w). (The effect of elastic supports in introducing changes in the 
moment is reflected in the value of the stability c nstant K), The change in (ie internal 
moment is evidenced by a change in the curvature (x m) and a possible change in the bending 
stiffness (EI). 


Referring to Figure 3,1. 2-1, we noie the equilibrium requirements. 


Mod 2 Mint = (1a) 


and 


-6(M_.) = 5(M 


ext int) 


evaluating the change in moments, we obtain 


~§ (Moet) = -6(Fw) = - [ P6w) +wi(6 ¥)| (ie) 
and a5 rao Rag 

S(M,u) = S(ED ay) = [EMO xy) + %y (5 ED | (td) 
Ae - | Fw) 1 wb r) | = | Hix) xa KD | (Le) 


_ Uf the load is monotonically increasing, then the column will fail when the buckling 
load (Fy) ig applied, No additional tncremertal load can be applied so that the buckling 
load is defined by setting 6F equal to zero, The above buckling criterion és equivalent 
to finding the load on a structure for which the deflection ts undefined (Figure 9, 1-1b of 
Reference 3-1), It ts also equivalsut to determining the load on the structure at which the 
ratio of 6F/6w is zero, where OF is the incremental load necessary to cause an incremental 
deflection (6w). The latter buckling criteria (@P) éw = 0) can also be employed to determine 
the stability of a column which exhibits creep, 


Setting OF = OQ and Fo = F, in equation (le) results in 


t 
- [ov + wd) | BLOX, * x45 EL (2a) 
ee ee OX ay oy bE] a is “s (1 4 SEI \ 
: 1 “" bw M éw jf 6y ng OX 5 i 
. OO 
ih Se et (=b) 


Thus the bucklime: load is a function of the ratio of the change of curvature to the 
change of lateral deflection and of the bending stiffness modified by an expression indi 
cating the cate of change of beading stiffness with cud cature, 


The first expression &x,¢, dw) is relatively simple te calculate if we assume that 


the deformation modes of the given structure do not change as the axial load {fs tnereased 
to the buckling load, Thus it is assumed that a pin ended column of constant bending stiff 
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ness, which buckles elastically in u sine wave, will buckle in a sine wave aven if it be- 
comes plastic, This condition is satisfied as long as the ratio of the bending stiffness 


distribution (oe) = Ela) EJI,) along the column does not change-significantly with the 


axial load. The stiffness can change due to plasticity but it is assumed that the ratio 
does not change significantly. This is consistent with our assumption that the bending 
stresses are small compared to the axial stresses so that the mn ‘erial moduli remain 
approximately equal (for constant cross--ection area). Thus the value of Oxy / dw is 


assumed equal to the value of the "linear column. This is the classical stability co- 
efficient found in various textbooks which can be expressed as the curvature to deflection 
ratio as indicated above or by the general expression of Eq. (1a) of Paragraph 2.3.1, 


2 
BX ay Fy 2 ? wy") dx 
= = = = Ay = Se = = “2. {3) 
bw EI ie (w,') dx L 


Values of K can be found in various texts (e.g. References 3-3, and 3-4), 


As an example for a pin-ended column of constant EI (i.e. p = 1) 


ot ain ME 
oa Se i eee ane eee 
vo TX a (4a) 
6 (sin) L L 


2 4} 2 
a Se * -5 == (4b) 
Gs cos 1) dx © 53 
t a 
Ken (4c) 
f _—_ 
: a i é6Et : ee 
The second expression | EI = EI + x,, ¢7— Jis evaluated by utilizing the defi- 
M 5 Xa 


nition of the bending stiffness derived in Appendix A of Reference 3-2 which approxir-ates 
the effects of the bending stress. 


From Eq. (Al0b) of Reference 3-2 


ren ets Ga) 
Els Esk a) (5: 
where x 
6 ee (ey oy) (Sb) 
fo So 
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Z 
and y= | E,pe/2E., . (5c) 


Substituting Eq. (8), (Ga), and (5b) in Eq, (2b) results in 


= K 2.2 

> Egohy (i - Gor’) (6) 

L 

The solution for the buckling load is now expressed in terms of geometry and boundary con- 
ditions (t, Lots K), the stress level at the neutral axis (Ego: En) and the stress distribution 


: x E 
(2 = —™ b- sty) Unfortunately the stress distribution parameter a is not speci- 
So 


fically defined but varies in a complex way with the applied load. It varies directly with the 
curvature (x ™ which increases nonlinearly with the load, inversely with the axial strain 


(€ é) which is defined by the stress-strain relationship (Figure 3. 2-1) and directly with the 


E 
stress distribution factor f ~ “* ”) which increase with load and curvature. Some ad- 
So 
ditional approximations are required to reduce the solution to a simple form expressible in 
terms of the initial conditions. Each of the approximations employed introduces small 
over- or under-estimates of the buckling load. It is hoped that the combined effect of all 
the approximations will esult in a reasonably accurate solution that considers the effects 


of eccentricity, 


&& 
FIGURE 3. 2-1 STRESS-STRAIN RELATIONSHIP 


The first approximation is to assume that the lateral deflection will increase as a 
roctangular hyperbola of the form 
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H} 


the lateral deflection 
fundamental deflection mode 


where w 


i] 


Wy 


ay = initial amplitude of the fundamental mode 


(The higher modes need not be considered in the buckling load problem since the 
column will buckle in the fundamental mode) 


Fy = buckling load for a linear or straight column 
K 
where EF = 9: Eglo (8) 


The secon approximation is to assume that y has the magnitude of unity in evalu- 
ating the value of a. Under the assumption of small bending stresses, the magnitude of 
y should be very close to but slightly tower than 1. 


The third approximation is to assume the average axial strain (€ ) to be equal to the 
value of F, 1/Eg of at buckling. 


The abeve three approximations are consistent with the assumpticn of relatively 
small bending strains and result in a slight underestimation of a and a resulting slight 
overestimate of Be 


The final approximation is to evaluate a with the stress distribution which exists 
at the cross section with the maximum lateral deflection. This assumption is made be- 
cause the stability is primarily a function of the square of the curvature (see Eq. (3)) which 
is largest in the vicinity of the maximum lateral deflection (for columns of constant cross 
sections), This approximation of the bending to axial stress ratio overestimates a and 
results in an underestimate of F, and compensates, to some degree, for the overestimates 
resulting from the small bending approximations. 


The approximations result in the following equation 
2/3 


E 
SAKE rs (“20) 10) 
oe 2 E Solo | 8\ a Egy / (9) 


These approximations cannot be rigorously justified but they result in an expression 
for the stability of the column that can be simply applied to the analysis of a very complex 
problem and which retains the correct sense for the effect of the parameters. A well con- 
trolled test program is recommended to evaluate the accuracy of the resulting formulation 
and to supply empirical correction factors, as needed. 


An examination of the experimental data of Reference 3-5 was conducted in Appendix 
C of Reference 3-2 and indicated good agreement with the results for specimens tested with 
eccentricities caused by eccentric loads or thermal gradients. ‘The predicted results were 
slightly unconservative. Better agreement can possibly be obtained between the experimental 
results and the theoretical predictions by statistically determining the best value for the ex~ 
pression v6 a,g/r). The graphs (Figures 3,3, 1-1 through -10) can still be employed by modi- 
fying the indicated value of ayo/r. Additional test data is necessary, however, to examine 
the reliability of the above approximations for greater eccentricities and larger buckling 
stresses than those reported in Reference 3-5 
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Using Eq. (9), the average buckling stress becomes 


pe ee ee ae oe rs (32) ( S22) ue) 10 
1A > 880. 42 NON y E. a 


So 


= ' ' 4 
- oO 2 a E 2/3 
[LE = = — as ve 9) ( 7 =) (ila) 
So L So 4 
0 F,/A Eo I 2 
Since €& = & — - = +. ee = K = {11b) 
So So 4 50 4 
€ a 7” 2/3 
a wie cz | 1 tas (“20) ( 3 sto )| | (llc) 
1 r So j | 
_ EA =O E 
Letting $= Gg ad @= ~ G 
fe) ) 
We obtain from Eqs. (9), (10), and (11c) 
F a 6 ‘ a.’ E,, \2/3 | 
ete an oat ge: ke cd ele ae Ae __To 
Ey Oo; ~ ey e - : Fe oy 2) | Gag) 


wr. tp) 


Su 
centricity and stress distribution, The effect is more pronounced with larger initial 
eccentricities or higher average stresses. It should be noted at this time that the factor 
was uviermined for sections which were symmetrical about the centroidal axis parallel to 
the bending axis (Qo = 0), Tne stability of the column should decrease slightly with in- 
creasing values of %0/Ig as indicated by a comparison of the bending stiffness in the 


plastic range for Qy) # 0 with Qy = 0 (see Eqs. (Ai@a) and (A10b) of Reference 3-2 | , 


The factor 1- 


2/3 
approximates the effect of the ec- 


3.3 ‘TECHNIQUE 


The final formulation for the approximate solution for the buckling of an eccontric 
column is presented by the non-linear Eqs. (9), (10), or (11) of Sub-section3. 2, [tis 
necessary to resort to graphical solutions in order to solve these equations. In order to 
simplify the analysis, a nondimensional solution is presented in Figure 3.3,1-1 ta -19, 
The graphs are based upon the stress-strain relationship presented in Section 3 of 
Reference 3-1, 
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3.3.1 Nondimensional Buckling Curves 


Equation (10) of Sub~section 3. 2 can be transformed to 


2) u (1) 
ae ) 
os Ee 7 ile 810 Baw 2/3 
E | Yr Ey 
A So 
PA Ree CARA 
se 4" o Ate] 
oA 0 
Utilizing the stress~strain relationship of Section 3 of Reference 3-1 7 (1 - B) re 
fe) 9) 


+B sinh °/o., , it 1s possible to calculate values of E/E, and En /E for various 


values of a/9, and given in values of 8. For a given value of the eccentricity ratio 
(219 /r) it is then possible to calculate the value of ©, which would result in a buckling 
stress ratio ©, /0,) , by direct substitution in Eq. (1). This relationship is plotted in 
Figures 3.3.1-1 to -10 for a large range of eccentricity ratios. The argument @ is a 
function of the material properties (Ey /o,) and the geometry and boundary conditions 
OKI At? ). The shape of the stress-strain relat‘onship is represented by the value of 8, 


The method of obtaining the material parameters from a uniaxial test is presented in 
Section 3 of Reference 3-1. The determination of the initial eccentricity which may be 
caused by initial waviness, eccentric loads, lateral loading, and thermal gradients is 
discussed in Paragraph 3, 3.2, 


The graphs for small eccentricity ratios should result in fairly good agreement 
with experimental data since the approximations employed should be satisfactory. The 
graphs for large eccentricity ratios are less accurate and are expected to be unconserva- 
tive, 
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BUCKLING OF ECCENTRIC COLUMNS 


FIGURE 3, 3,1-2 
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FIGURE 3.3, 
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3.3.2 Initial Eccentricity 
The initial eccentric.:y must be evaluated in order to use the analysis curves. 
Paragraph 2.3.3 indicated that the amplitude of the fundamental modes 
can be obtained by a weighted integration of the actual lateral deflection or by matching 
the deformation at discrete points. 
Assuming that it is possible to express the initial eccentricity (wo due to shape, 
axial load location, lateral loads and temperatures) as an analytical expression, it is 


possible to determine the magnitude of the first mode by utilizing the orthogon: ty be- 
tween the lateral deflection and curvature (Eq. (2) of Paragraph 2.2.1), 


Thus, if 


Ww. (la) 


then from Eq. (2) of Paragraph 2. 2.1 we obtain for constant beiding stiffness (c=1) 


4 
i wx, dx 
P o 4 
Bitte Sager (1b) 
{w, », dx 
1} 1 
0 
The following formulae are applicable to columns of constant bending stiffnes .° 
when j 
wy = 5 mig (2) 


then for pin ends 


a 


10 7 20m, Sq, L) (3a) 


while for clamped ends 


ay = 2Em, CG, 1) (3b) 
when 
1 
My = op Em e (4) 
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then for pin ends (Reference Eq. (a) of Paragraph 2. 3. 3.1.1) 


Noo > 2 apeyasey <: ° <8) 


and 


a 


m 


2 22 Garten [sa+2, 1 -S(,n) | 


while for clamped ends (Reference Eq. (4c) of Paragraph 2.3.3.1.1) 


m, : 
Wy = Loa (OP -1 8 - ee) 
and 
a = 22g [cow H-508, 9-a-po@, 0 | 


Values of S(j,1) and C(j,1) are found in Tables 2,3.3,1.2-1 and ~2 and are 


summarized in Table 3, 3. 2-1 below. 


TABLE 3,3.2-1, 


. 95066 
. 07599 
, 08592 
. 08815 
. 08669 
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FOURIER COEFFICIENTS FOR POLYNOMIALS 


(5b) 


(8b) 
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For a pinned end column of constant EI and constant linear thermal gradient 


0 
7 m_&é 
“a = shor = ~ 2 we obtain from Eqs.(Ga) and (6a). 
c h rh 
m m 
te) 0+2 O 442 
Yo = (rom ean a ae 
m m 2 
: Op ge - o _t' AaT 
ve woth = =~. 68°55) = += ee 
4 rc ‘ . 1 1 | 
Bee 419 = 2 ae [8@%)- SG) | 


m a 2 
. r) a _ ,129L° AaT 
ore ap = 2 = is 1893 - 3183 | = -,129m) = - oo (7) 


it 


A second method of obtaining the initial eccentricity is by matching the displace- 
ments at a discrete number of points, This requires the sotution of the set of simul- 
taneous equations w{é) = Da, w, (€) for the value of a5 and is described in Paragraph 


2,3,3.2. The approximation a,, = We (1/2) is a solution where only one point, the de- 


10 
flection at the center of the column, was matched. The accuracy would increase with the 
number of displacements which ace matched, although the matching of the center deflec- 
tion of a pin ended column with a uniform lateral load or thermal gradient results ina 
satisfactory determination of the initial eocentricity. The accuracy also dep ds upon 
the form of the lateral deflection. For example, a matching of the mid-lengt. s deflec- 
tion for a pin-ended column resulis in a more accurate amplitude of the fundamental 
mode caused by a uniform thermal gradient (parabolic) than by a constant eccentricity, 


3.4 SPECiAL CASES 
There exist special situations for which the exact solution is known or for which 
engineering approximations have been accepted because of experimental data, The 
plausibility of the approximate formulation, presented in this section, is reviewed by 
degenerating {t to those situations with which the analyst has had some experience. 
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3.4.1 Linear Material 


E 
A material whose modulus is independent of the stress levels e =0, E= Ee 


a Eg) is described as a linear material. The value of (1~ E Eq) is identically zero 


To 
and no reduction of bending stiffness occurs, The eccentricity (a, of r) does not affect the 
stability and buckling occurs when the average axial strain attains a value of K (=) (with 


the load equal to Fe = * El) provided no fiber is stressed to its ultimate strength. 
t 


This corresponds to the classical "Euler Column" whose buckling load and strain are 

well defined because of the linearity of the material beyond the buckling stress. A column 
with a large slenderness ratio (t/r) would behave in the manner described, Although the 
nondimensional stress-strain relationship includes the linear case (letting B = 0), it is 
recommended that the actual value of 8 for the structural material be used 2ven for the 
case of low buckling stresses, This will permit an approximation of the effects of eccen- 
tricities in causing extreme fiber stresses that may be beyond the "proportional limit" 
even when the average stress is very low, The deviation cf the stability of a structure 
composed of material which is non-lin ar to a slight degree from that of a linear material 
will be insignificant for small eccentricities but may become significant for large eccen- 
tricities, 


3.4.2 Perfectly ‘ caight Columns 


The perfectly straight column does not bend when subjected to an axial load, This 
condition is virtually impossible to attain experimentally but has some theoretical value 
as a mathematical model. The column will be stable is long as the axial strain remains 


below the critical strain K — since 249 is zero, Any lateral excitation of the column 
e 


below this value of the critical strain wiil dampen out whereas it will magnify and become 
excessively large when the critical strain is reached or exceeded by the column. 


3.4.3 Small Initial Eccentricities 


The effect of the eccentricity upon the stability of the column is small when ie 
eccentricity ratio (a,,/r) is small. An upper bound upon the stability load will always 
be F; which corresponds to the employment of Eggo as the effective buckling stiffness 
of the cross section and a critical strain of K r2/t2, 


Under the conditions of continuous loadings (Eo = Many) and small eccentricities 
(2)? “ae Py oe: ee em 
(LS b= ), the lower bound of the stability load will be F,{ > => EL. 1) which 
fo # 1 Bog Fig To 0 
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corresponds to the employment of Eno i, as the effective buckling stiffness of the cross 


Ex, y* \ 
section and a critical strain of | =—~ | \K “> 
Eso vu 


This lower bound is satisfactory provided the eccentricity is not too large. The 
bending stresses become more significant in ueterraining the bending stiffness as the ec- 
centricity increases, An eccentricity ratio exists for each column beyond which the tan- 
gent modulus stability ceases to be a lower bound. This eccentricity ratio can be obtained 
by equating the tangent modulus load to the eccentric column load, 


a 1-E E 
Equality results when ~1o = = 
E E 
A value of (a, /r less than 1. “to/*s0 results in a tangent modulus load 
that is conservative with respect to the eccentric column load (i.e., Enolo s ET) whereas 


a larger eccentricity ratio will make the tangent modulus load unconservative, The value 
of the critical eccentricity ratio depends upon the material and the slenderness ratio of 

the column, This is illustrated in Figure 3.4.3-1 which indicates the effect of eccentricity 
on the stability of the columns described in the illustrative problems, The tangent modulus 
loads stops being conservative for eccentrivity ratios between , 2 and .5 for the columns 
analyzed, and the critical eccentricity ratios are indicated in Figure 3,4. 3-1, 


3.5 ILLUSTRATIVE PROBLEMS 
The computational techniques are illustrated in the following problems; 


Find the buckling loads of the following pin ended columns of constant cross section 
jilustrateu in Figure 3.5-1. 


ry ci a) 4 = 100" (Slender column) 
SEE = oa " ig c 
T, = 34°F ep ; FS aeTeee b) 4 100/.,/3 =62.1" (Intermediate column) 
04 ee = ce) 4 = 100/10 =31.6" (Short column) 


x Constrained to bend 
about this axis 


= 


SEC A~A 
FIGURE 3.5-1 COLUMN SUBJECTED TO LOAD AND TEMPERATURE 
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The temperature is assumed constant on each face of the column with a linear 
gradient through the thickness. Material parameters are assumed for 2024-T86 at the 
mean temperature of 300°F to illustrate the computational technique. The actual values 
must be determined from standard material tests. 


Material Properties 
The following material properties are assumed: 
12(10) in/in/*F 


(¥ = 
Tes a 
EA = 10 psi 
3, = 5000 psi 
Bodo: 


Initial Geometry 


‘Ae OVE 0,04 999+ 10) = 9 in” 


ae 2 3 4 
1, = 2612)" + Fp ¢0G.9)" = 2.30 in 


ro 
r? = T/A = 2,2 in” 


yr = 1,15 in 
Fa xr Fa . x*? 10) 22.2 4.4 
Re SG ee ee ate ea ge ote 
o AL 0 2? 5(10) £2 (4/100) 
2 
q= KS = ee ee 
a L 4 


Initial Eccentricity 


Because the temperature distribution is linear through the thickness we can utilize 
Eq. (4) of Paragraph 3. 3.2. 


re amen 
eh ae ea gs ee " 
10 h h 
r 
819) 129 w (T,-T) 0” ee , 
mE SN Teh = Az0cnsy toy” (948-257) _ og( 4 
_— rh = SIiyaay OD 
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From standard reference texts the lateral deflection due to a uniform load q is 


4 
Wy = gear e- 26 + 84 


From Eq. (3a) of Paragraph 3, 3.2 


af 
Qi) g 22 HiT [ sa,1) ~ 28(2, 1) + S(4, | 


From Table 3.3, 2-1 


4 4 
= gt... as | — 1 qt 
@1q)g = 2 Say [- 818. - 2(. 1248) + (08814) | = .013066 - (2) 
(a, ,) 4 4 
10g _ ee = .08 (;*.) 
i (10) (1. 8)(1. 15) 
It should be noted that 
5 as Lge qt! 
7 t = % =e = be MERE ey ass 
Wg (at x = 4/2) = yoy 4 01302 “rym. 013066 ET- = (Ay 9)g 
and 
ae =) gage OE G2 AaT ,2 ; 
(W dp (at x = L/2)= .125 h 4° x .129 ae = (@)9)p 


This indicates the magnitude of the errors which can be introduced by approximating 
the amplitude of the fundamental mode by the deflection at the center of the column are quite 
small for a pin ended column subjected to uniform lateral load or thermal gradient, 


The initial eccentricity ratio is then calculated as follows: 


ia) 5 a ia sagt el 


r r r 100 i100 
For 2 
} 

4 = 100 : -#f 28 + .08 - 236 
a 16 

4 = wo//3 ; FP OS 4 OUR. ee TUL 
aa Qo 

L -- 100//10 : “Eo ORB 0008 = , 0288 
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Stability 
_ Referring to the appropriate graph of Figures 3.3.1 results in the determination 
of (,/0,) for the given velues of 8, ©, and (24 9/,) The effect of the eccentricity, which 
is expressed in Eq. (11d) of Paragraph 3.2, can be evaluated by comparing the given ralue 


of € to the value 6 which corresponds to 6,4, } for a straight column (Figure 3.3, 17-1 


E 


with a = 0). This value of = #3 €, , corresponds to a material and geometty 
10/r a, 2 


parameter for a straight column which would buckle at the same stress as the eccentric 
column and is a measure of the average axial strain in the column. 


a) For 4 =100, E, = 10" g, = 5000, and B= .0001 we obtain 
419 
@= 4.4, €, =. 0022, and — = 36. 


From Figures 3.3. 1-7 and ~8 


G 

1 

Gre 4.4 

16) 

an = 4,4(5000) = 22000 psi 

F = G, A= 22000 (.59) = 130004 


From Figure 3.3, 1-1 
b= 4.4 


@/b = 4.4/4.4 ~ 1 (no effect of eccentricity because of the 
large slenderness ratio} 
and x 
C2 «, (ae) = .0022 
b) Similarly for 4 = 100//3_ 
= 
@ = 13.2, €, = 0066, and —— = .10l 


From Figu.e 3,3. 1-6 


ec] 


S, = 9.75000) = 48506 psi 


es) 
It 


48500(.59) = 285004 
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and from Figure 3, 1. 1-1 


$= 10.5 $/@ = “2 = (80 


J. €, = ,0066(.80) = ,0053 


ce) For’ = 100//10 


ia) 
@ = 44,  €, = .022, and —— = ,0288 


oF % = 43,4. 


Cs = 1£,1(5000) = 65509 psi 
Fy = 65500(,59) = 28500# 


The eccentricity ratios Qosp were sufficiently close to available values that there 
waa no need to interpolate between graphs. Plots can be made to show the variation in the 
stability of a given column with the initia! eccentricity awl can be utilized if interpolation is 
required, These plots are shovn in Figure 3.4.3-1 for the three gong, intermediate and 
short) columns ana'yzed above. The values corresponding ie an effectiva rodulus of Bs 


: » Qe R : : 
sre obtained from Figure 3.3. 1-1 ee = 0) and the values corregponding to an effective 


: 
modulus of Ep are obtained from Figure 9,2.1-2 ut Reference 32-1. The pote indicate 
that El = EL. Js always uneconservative and EL = EL. 7 becomes unconservative in 
Soo : Yoo 
Xp 
the vieinity of .2 > ae > 9 for the analyzed columis. ft should be noted that the re- 


ductfon in the stability of the long column is least white the reduction in the stabilfty et the 
intermediate column is greatest with the shert column dfected to an intermeduite deyres. 
This is because the reduction in the stability of the eceentric columm depends both upon the 
eceentricity and the stress levels attained by the columal The long colunin bas the largest 
eecentricity ratio bul has very low axial stresses, the short column has very hich axial 
stresses but very small eceentvicities; while the intermediate column has high stresses 
and moderately high eceuntricities, 


It is interesting to nefe the contributions of the various factors in reducing: the 
fe K m1 \ 
ii Ee Ta as Vhe factors redueing the 
LE CO. 

stability stiffmess ean! rourhly divided into (wo parts. The tirst part represents the re: 
duction in the axial stiffness due to the plasticity caused by averave stresses which are above 
the proportional limit of the materiai. The second part represents additional reductions in 


stability of the column below the “Euler Loud" 
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the stability stiffness because of the eccentricity causing a shift in the neutral axis and 
a rate of change of the bending stiffness. 


Equation (2b} of Paragraph 3.2 can be transformed to the following approximate 
form in order to evaluate quantitatively the destabilizing effects of eccentricity and plas- 


ticity. 


#y 
Pe 


where 


rs 
ES. I 


EI | x aEl Esl BY) (1 x El 
EX, Eg, Eool, 3% \ EI, Bool, Fol, OX 


ig the buckiing load of the eccentric column 


is the Euler buckling load 


ig the factor representing the reduction ia axial stiffness 
is the effect of the shifting of the neutral axis 


is the eifect due to the changes in the bending stitfness 
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TABLE 3.5-1 STABILITY RATIOS 


“LENGTH OF COLUMN 


| 4=100" | t=100/,/3 | 4=100//10 | 


DESTABILIZING Pi. NOMENON 


Reduction in Axial Stiffness 
(Plasticity Effect} 


E, I F 
soo 1 
El * F. 1,00 .73 _ 30 
° E 


Further Reductions in Stability Stiffness 
(Eccentricity Effects) 
a) Shift of Neutral Axis 


as ee ee 1,00 (93 955 
m1 
So°o 


b) Rate of Change of Bending Stiffness 


1- -% Met 1-4@7 Pr" 1, 00 87 91 
S00 e 


c) Critical Strain Ratio- Total Eccentricity 
Effect | 


2 


2 1- 6o'r? = 1.00 | — .80 865 


Stability Stiffness Reduction 


F FL, OF 
aE he 1,00 .o8 . 26 


1 
Pa 
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SECTION 4 - AXISYMMETRIC LARGE DEFLECTIONS OF 


CIRCULAR PLATES SUBJECTED TO THERMAL AND MECHANICAL LOADS 


4.1 SUMMARY 


This report is concerned with the nonlinear axisymmetric analysis of circular plates 
with in-plane edge restraint Both temperature and mechanical loads are accomodated as 
an extension of investigations performed for the isothermal mechanical loading problem 
(References 4-1 through 4-4), An exact mathematical formulation within the framework of 
the v. Karman large strain-displacement relations (Reference 4-5) is developed, The 
equilibrium equations and boundary conditions are then derived by utilizing the calculus of 
variations for arbitrary axisymme'rical temperatures and normal distributed loading. The 
satisfaction of equilibrium and compatibility equations requires the solution of two simul- 
taneous nonlinear ordinary differential equations subject to the prescribed houndary condi- 
tions. Analytical solutions of such equations are apparently not possible and therefore num- 
erical procedures must be employed. 


A finite difference procedure utilizing "relaxed iterations," developed by H, Keller 
and E, Reiss (Reference 4-4), and employed by them for the solution of isothermal problems 
with apparently unlimited load parameter ranges, is used here fur combined thermo-mechani- 
cal problems. Numerical results are presented for the special case of a simply supported 
circular plate with radially immovable boundaries, subject to a uniform pressure and an 
arbitrary temperature variation through the thickness (no planform variction), These re-- 
sults have been obtained for a large range cf temperature and load parameters, However, 
because of space limitations, only a limited amount of data is presented in this report. 


4,2 INTRODUCTION 


One of the basic assumptions of the classical linear theory of plates ig that bending 
action does not induce significant midplane stretching, It is further assumed that stregses 
and deformations produced by loads and restraints in the midplane are superposible cn the 
bending solution, Thus, coupling between the two offects is not accommodated by the clas- 
sical theory. When the deflections are not small compared to the plate thickness, midplane 
stretching is no longer insignificant, resulting in a nonlinear interaction between bending and 
membrane stresses, Therefore, large deflection theory must be employed. 


It is the purpose of this report to Investigate the axisymmetric large deflection prob- 
lem for circular plates, 


The general formulation presented considers arbitrary axisymmetric temperature 
and pressure variation, where the von Karman large strain-displacement relations (Refer- 
ence 4-5) are utilized, These assume infinitesimal strains and finite but small normal de~ 


2 
flections ( eu is of the order of the strains, but small compared to unity), The remain- 


ing assumptions are those of classical plate theory. This formulation is more complete 
than the conventional linear theory in that the results are valid for deflection magnitudes 
several times the plate thickness, Moreover, buckling and poatbuckling behavior are om- 
bedied in the analysis. 
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4,2 (Cont'd) 

Numerical results in nondimensional tabular and graphical form are given for a 
simply supported circular plate, with full boundary restraint to radial movement, sub- 
jected to uniform pressure and arbitrary temperature variation through the thickness. 


The following symbols are used throughout this section: 


b Plate radius 
h Plate thickness 
r Radial coordinate 
r rs nondimensional radial coordinate 
u Midplane radial displacement 
uw, Midplane radial displacement at plate edge 
uw Radial displacement 
w Normal deflection 
W M4 , bondimensional normal deflection 
Zz Thickness coordinate 

Eh? 
D aE flexural rigidity 

12(i-v") 

E Young's medulus 
K Relaxation parameter 
M_, M. Radial and tangential bending moments, respectively 
M. Mm Nondimensional radial and tangential bending moments, resp.rtivcly 

h/Z 
My | RaTzdz 

-h/2 
My. Nondimensional form of M.. 
Nw N, Radial and tangential membrane forces, respectively 
NWN, Nondimensional radial and tangential membrane forces, respectively 
4,3 
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h/2 

Ne J EaTdz 

-h/2 
Ne Nondimensional form of Neo 
q Normal pressure 
Q Nondimensional normal pressure 
T Local temperature with respect to an unstressed and undeflected datum 
U, Strain energy density 
Vv Total potential energy 
a Coefficient of linear thermal expansion 
B Slope 
Cre, Radial and icngential strains, respectively 
en € : Midplane radial and tangential strains, respectively 
r Elastic in-plane edge restraint 
v Poisson's ratio 
¢ ®, nondimensional stress function 
o; Finite difference value of ¢ at the i'th grid point 
¥ Stress function 
Oy,» Radial and tangential stresses, respectivaly 

D Igquy? 
6 i joa v) 6 
8; Finite difference value of ¢ at the i'th grid point 
4,4 


4.3 BASIC EQUATIONS 


4.3.1 Sress-Strain - Displacement Relations 


The von Karman large strain - displacement relations for the axisymmetric strains at 
any point in the plate are given by 


2 
_ due, 1 (dw 
Ss * ar * 3 ($ 

. (1) 


yt 


e 


where u*(r,z) is the radial displacement of the point and w(r) is the deflection normal to 
the undefiected midplane (Figure 4.3. 1-1). 


FIGURE 4,3.1-1 GENERAL AXISYMMETRIC DEFLECTION OF A TYPICAL POINT IN 
THE PLATE 


Assuming plane sections tc remain plane and normal to the deflected middle surface, 


dw 
u* =u~2 se (2) 
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where u=u(r) is the radial displacement of the midplane. From (1) and (2), 


Jai, & diy 28 dw? 
&. ~ ar eae. ee 
(3) 
-u_% dw 
%& “yy 7 r dr 


The last (nonlinear) term in the first of (3) dos not appear in the conventional small 
deflection theory. From Hooke's law (neglecting normal stresses in the thickness direction), 
including the temperature "T'', there results 


E 
as 5) [e+ ve, - (1+) at | 
1-v 
(4) 
ET 
ee” [ep t+ ve, - Grr) | 
-Y 
An integration of (4) through the plate thickness yields’ 
h/2 r 
_ Eh re) o_ (ity) 
N. = | 0, dz - ae e. +6, Eh Nip 
-h/2 
(da) 
h/2 
_ Eh |oo,,,0_ (Ly) 
Nf ode = Sr ia ines ge 
-h/2 _ 
where the midplane strains are given by 
o, du, 1 (a : 
op de dr 
(4b) 
Bi ct ee 
"to eM 
t ia what follows, it is assumed that aT =aT{r,z) but that E (and hence D) is constant. 


roe 
toy) 
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and 
4/2 
Nip = [ EaTdz . (cj 
-h/s 
Multiplying (4) by z and then integrating through the thickness yields 
h/2 
2 M 
z --p|/a¥,% dq, Tt _ 
aes Cre ar" +r a” D(i-v) 
~h/2 
(4d) 
h/2 
2 M 
a -~p}i dw dw qT 
M, ae 0, 2 dz = -D - apo a + Di-v) 
-h/2 
where 
h/2 
Mp = EaT 2dz (4e) 
~h/2 
and 
gh® 
DSS re 
12(1~v") 


4,3,2 Strain and Potential Energies 


The strain enorgy per unit of volume is given by (Reference 4-6): 


1 P 
U, = 3 Exe + & O gee Ue (Pa 0) | 


Substitution of the stess-strain relations (4) of Paragraph 4.3.1 into (1) yields 


(2) 


eee 


E 2 Di sel 5 Fahy a2 
a aes [ + & + 20 6.8 - 2(1+v) (e+ 6+) aT+ 2(1+v){aT) | 


For rotationally workless restraints at the outer boundary, the total potential of the plate 
and external loading system is given by the following equation: 


A a teene eee e ae 
2 2. SSCP TT Setters, es 
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b  h/2 b 
2 
= an { r Uj dzdr ~ an { qwrdr+ mb ru, (3) 
o 8 6«(-h/2 re) 


where q is the normal pressure, \ is an elastic restraint per unit circumferential length 
to radial displacement of the boundary, and uy, = ul che From Eqs, (3) of Paragraph 4, 3.1, 


and (2) and (3) the potential energy in terms of the displacement components and temperature 


becomes 

* - 2 

v= [eli [B25] cy a [e 1G) 
) 


dw\? fid y av dw déw\ bh? 

— + \ — me + ——_— —_— —_— 

dr? xr dr r dr ar” 1 

(4) 
h/2 
M. 2 

pe dz+ | dw, i dw 
+ 2(1+v) (aT) dz + 7 : 5 tt | 
-h/2 -! 


ey, 
mi 


2 vs 
Tldu, 1 fdwy .u 2 
me ape oe ge dr - 2n [qwrdr + mba (uy) 


4.3.3 Governing Differential Equations and soundary Conditions 


The equilibriwn equations and "natural" boundary conditions are now obtained by 
making the potential energy stationary with respect to variations of the displacements w 
and u; i.e. bv ES gv ~Q. The firsi of the variations (6,V = 0) yields the equilibrium 
equation 


4 Ld (i4 
DY Wes ae Ng) = q~ l-v (2) 


and the following boundary conditions: 


Or Sa0 
Assumed regularity of the siution at the center requires that M_ and M, be 
finite. This implies that eer and M =M,. Ge) 
dr r t 
4.8 


4.3.3 (Cont'd) 


(2) r=b 
Gite eg dw 
(i) w prescribed or - ar [pv Ww +t |e N. ae 0 
(2b) 
2 M 
ii ' GoW Uh ke 
(ii) w prescribed or D . 5 ae ie hg 2 
The second of the variations (6,4 = 0) results in 
N, - a rN )=0 (3) 
t~ dr @N,)= 
with boundary conditions: 
(4) r=0 
N.. and N finite. This implies that u=0 (hence N.. = N,) : (4a) 
(2) r=b 
N.. 
ut) =0 or u prescribed. (4b) 
The two equations (1) and (3) contain three unlmown functions of r; i.e., wy, 
N. and N.. A third equation is obtained from the necessary condition that these three 
quantities yield a set of single-valued displacements, u and w. A statement of this 
requirement ig obtained by eliminating u from Equation (4b) of Paragraph 4.3.1, which 
results in 
0 ° fe) 
ae ec Bs 
dr r oe ar \dr/ * (Sa 
Substituting (4a) of Paragraph 4.3.1 Into (Sa), 
dA ncn OD ong EE YY a (3b 
ap he ON Ny) pr NN Y oy dry ae 


Equation (3) is automatically satisfied by a stress function 4 defined through the 
relations 


y= tN. 

(6) 
dy _ 
dr Ne 
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Substitution of (6) into Eqs. (1) and (b) results in the following coupled set of nun- 
linear differential equations 


1 
pvw-2 2 # - Ey VM, 7) 
and 
2 wn? 
r dr = dr ~T 
dr 
where 
ee ea 
ee r dr * dr 
e 2uy < 
A more convenient form is obtained by introducing the slope § = ay , which yields 
aM 
d d id ; iz oe oe “T 
Dar tar r de (B)> ap 8) = a Ty a( dr @) 
and 
a’y yoy, Eh 2 eae 103 
dee ob EE or FP ar Ny oy) 


Equation (9) may now be integrated with respect to r, resulting in 


r 


ae 
GP DGB Pod Be od if eee See 
Dl (= &) ok ae | eee Mae ae Mee (11) 
QO 


The Eqs. (10) and (11) are to ba solved subject to the boundary conditions (obtained 
from (2) and (4)) 


(1) Bo) = 
(ii) 8(b) prescribed ae si rpe) (12) 
or : 
da v + L = 2 2 
Oe gee ae ae = (no radial bending moment) 


(i) yfo) = 0 (since NS ig finite at r = 0) 
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(iv) a + ( Bhs g 2} yt Ny = 0 (elastic radial restraint "\" at r=b ) 
e bo) / mes 
[ st we V+N, | prescribed (known radial displacement at r=b ) 
| dr r Po cb 


The above differential equations and boundary conditions can be expressed in non~ 
dimensional form as 


6 8 94 “yd 
Cie = oar ea =rQ(r) M, 

J r r 

13 

i a 4 (13) 
ht ¢ a OL tee a N ~ 

pee : a 

r Yr Yr 

where: 

peo ie pel 

b ? 
(+= a 

dr 7 
tb 
Ag 
a - h Y bil-v) 8 
rm Shae a 2, = 
Qt) s FB bectee"y; aE oy (@) Ede (14) 
Sa Ae a a 
Eh 1 
oO 
et So) 
ae Dh 1-y 
a 

Nb 
pe ie 
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The boundary conditions become: 


(1) e(0) = 0 
(2) 9(0) = 0 
(3) 9(1) prescribed 
@(1) + vO{1) + Mp (1) = 0 4a 


(4) acy +( FP -v) o(1)+N,, () = 0 


or 
(1) ~ v g(t) +N..(1) prescribed 


The set of nonlinear differential Equations (13) and accompanying boundary conditions 
(14a) are not amenable to analytic solution. Numerical procedures must be used to obtain 
solutions for specified values of the parameters. A finite-difference approach is presented 
below together with numerical results. 
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4.4 NUMERICAL INVESTIGATION 


4,.4,1 Finite Difference Procedure 


The procedure developed in Reference 4—4 is employed here f the combined thermo~ 
mechanical problem, In particular, we consider tne case of radially immovable edges 
(\ > © in (14a) of Paragraph 4.3.3) with simple supports for berding where q, Mis and N,, are 
constant (uniform pressure and temperature which varies only through the thickness), For 
this problem, (13) and (14a) of Paragraph 4.3.3 reduce to 


Ld =09+QF" 
g (1) 
La-=-9 
and 
8 (0) =0 
§ (0) = 0 
6 (1) + ve()+M,, = 9 () 
T 
@ (1) - vd (1) + Ny. = 0 
where 
eee 
dr r dr 


4 3/2 
oR) [aw] 


and the cther nondimensional quantities are as defined previously. 


A central difference representation of the differential equations (1) vields the following 
foi the interior points: 


(td te inet} (3) 


where for an arbitrary function ¢, 


D 
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Bl ee Sei | ke ; im : : 
Eds ete Leena. 714 | ~ poe LIG GG | (3a) 


8 
MN 


number of subdivisions 


~ at aa: 
a or ee 
m m m 


N 


The finite difference form for the representation of the boundary conditions may be written as 


84 = 0 
#5 = 0 : 
b 
mé = cas j6(1 +") 
a m-l Dh Nin Topas (4) 
m (m+ Vv) 
N..b" 
‘ 7 
[mo = ae 
eo =} m-1] D 
m L m-v 


where backward differences are used for the end point r= 1, 
Equations (3) and (4) constitute a set of 2(m+t) equations in the 2(m'1) unknowns 
8, . (i° O,1,2...m). However, these algebraic equations are both coupled and nonlinear, 


and cannot be solved in ciosed form. A "rele :.4 iteration" technique is employed, in which 
the iterative forms of equations (3) are written as 


ian ° [3] ne esi (PP) 
le] | x(a] [2 : « }[e,] (5c) 


Giz 1,2,3...m-1) 
where K is a "relaxation™ parameter. 


A typical iteration, starting with the n'th set of iterates 18, | is as follows: 
Po) : 
(1) Substitution of 1 8 into (Sa. and making use of the second and fourth of 


Eq. (4) yiet ds a set’ of {m-1) tri-diagonai, linear algebraic equations from 


which the quantities | ¢, | are determined. 
nel 
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(2) Substituting fe, |» [os Ia into (Sb) again yields a set of linear equations from 


feed 


* 


tree F 
eel for the next set of iterates are determined. 


which the provisional values [6 


(8) The actual value of the next set of iterates cae is obtained from Eq. (8c) 


where a suitable value of "relaxation" parameter is employed. 


(4) _Iterations to convergence are performed, 


4,4,2 Numerical Reyuits 


Based on the procedure indicated, results are presented iu both tabular and graphical 
form (Table 4.4,2-1 and Figures 4.4, 2-1 through 4,4, 2-6) for the simply supported solid plate 
with radially immovable edges. 


The nondimensional deflections (w), membrane forces (N.. and Ni), and bending moments 


(M., and M,) presented in the graphs and tables are defined as follows: 


Wine 
h 
Nb? 
N - r 
r D 
_ Nb 
N, Fg ee ee 
2 
Mb 
MF - yr »f 6(1-v ) 
(4) 
M,b® pecs 
M, “ “Dh jf 8G -v } 


and as defined previously, 


Pe. he 
5 
9 ote h/2 
zr OCOD (itv) - 
Mp pr ipa? feted 
+h/2 
5 h/2 
bcs Fb , 
ie mer Oe) ne? 
Nv. p | aT. 
+h/2 


4.15 


4.4.2 (Cont'd) 


Stresses (in nondimensional form) may be obtained from these formulas {Reference 4~7): 


bho, _ iag f Mp ~ M, a Ny ore bp? ~ 
D ae erence ia (i+v) —y oT + Ny 
= = 
[6a-») 
(2) 

pho M,.. - M, N 2 
t _ 122 7 t 7 b — 
a Se ey Ee 


/ 7 2 
j6(-v") 
A discussion of the numerical results follows. 


(1) Thermal Buckling 
The thermal buckling problem of a circular plate due to an average elevated tempera- 
ture through the thickness (proportional to N.,), where M.,, ~ Q = 0 and radial edge displacement 


is prevented, is equivalent to the mechanical buckling problem of a plate having an edge thrust 
corresponding to a prescribed edge displacement, It Is shown in Reference 4-2 that buckling 
can occur only when the edge thrust exceeds the lowest eigenvalue of the linearized buckling 
problem, The critical 'hrust corresponding to this eigenvalue ig given by (Reference 4-5): 


2 
(ae. 2.05) D 
-[-,; ° See 3 
cr beh 
Since, for the thermal problem, the thrust (up to buckling) is given by 
o ee aE r (4) 
r h(l-v) 


then from (3), (4), and the definition of Ni we find that 
N.. | 2, 94 
L Vor 


The postbuckling behavior of the plate is given by the nonlinear analysis; the numerical resulta 
are presented in Figure 4. 4, 2-1 and the first five sub~tables of Table 4.4.2-1. The table eni- 
ploys a floating decimal number system which is ‘o be interpreted ag shown by the following 
examples: 


O.5159E 00 (0.5159) x (10°) - 0.5159 
0.5159E 02 (0.5159) x (10%) 51,59 
1 
) 


0.5159E-01. (0.5159) x (10--) > 0.05159 
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(b) RADIAL MEMBRANE FORCES 


(a) DEFLECTION 


“120 ee ee: 
0 2 4 6 8B 10 


7 t 


(c) TANGENTIAL MEMBRANE FORCES {dq} RADIAL BENDING MOMENT 


FIGURE 4.4,2-1 NONDIMENSIONAL THERMAL BUCKLING AND POST BUCKLING 
BEHAVIOR (M,, = Q* 0). 
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Figure 4.4, 2-1 shows the variation of the nondimensional deflection, membrane stress 
resultants and bending moments with the nondimensional radial coordinate r, as Na varies 


from 0 to 100. It is noted that for 0< Ny < 2,94 the plate is undeflected, and the two- 


dimensional linear elastic solution holds, yielding stresses which are conctant over the plan- 
form and zero bending moments, For higher values of Nes Figure 4,4,2-1a shows nonzero 


deflections that increase monotonically with increasing N, For low values of No the 
quantities N. and N, are compressive (Figures 4,4, 2-1b and c) throughout the plate. With 
higher values of No (25, 50, and 100, for example} deflections in the central region appear 


to be restrained by tensile membrane stresses, while in the vicinity of the plate edge the mem- 
brane stresses become compressive, Figures 4.4, 2-1d and e indicate that for small N.p (but 


above that causing initial buckling), the maximum moments in the radial and tangential directions 
occur at the center ofthe plate. With larger Ro maximum moments are awav from the center and 


approach the outer edge as No increases. The radial moment, in particular, exhibits astrong 


"boundary layer" effect (sharp gradients near the edge) in meeting the condition | M, [| =, 
r=] 


(2) Bending Due To Combined Thermal And Mechanical Effects 


Starting with the sixth sub-table of Table 4.4, 2-1, nondimensional numericel re- 
sults are listed for the thermal bending problem (Q = 0), corresponding to the following com- 
binations of thermal parameters: 


M.. = 60, 100, 1000, 2000 


No = 0, 2.94, 5, 10, 25, 100 
For a qualitative discussion of the results, we refer to the graphs of Figures 4,4, 2-2 
through 4,4, 2-6 in which mechanical loads (Q 4 0) as well as temperature are considered, 


(a)  Deflections (igure 4. 4, 2-2) 


For low pressures (quantities proportional to Q) the deflection in the plate interi: r 
becomes constant for large temperature differences (quantities proportional to My). This 


flat deflected shape is maintained almost to the plate edges where the bounc .ry requirement 

of zero deflection causes a boundary layer effect {n which the deflections drop precipitously, 
Thia effect becomes less pronounced as the pressure loading increases. However, it is 
interesting to note that even for high pressure loading (Figure 4,4. 2-2d), the plate still tends 

to flatten out in the central region as the temperature difference increases and, contrary to 
what would be expected, the deflection at the center does not in all cases increase monotonically 
with the temperature difference. ~ 


{b) Membrane Stresses (Figu 


With increasing temperature difference between the plate faces, the radial and 
tangential stresses tend to become conatant and equal to each other in the plate interior, 


reg 4.4, 2-3 and 4.4, 2-4 


(a) 


{c) 


FIGURE 4.4,2-2 NONDIMENSIONAL DEFLECTIONS DUE ‘TO TEMPERATURE WITH 
AND WITHOUT PRESSURE LOADING 
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FIGURE 4.4,2-3 NONDIMENSIONA!] RADIAL MEMBRANE FORCES DUE TO T™MPERA TURE 


xil 
+t 
a 
A 
._o 
to) 
(o} 


SE ae a ciets Coane 
i 


BIGURE 4,4, 2-4 


560}-— Bal 
| | i ae : 
s | | | “tows 7 ] 
480/- | ewe 1 
[Sic Sete = 
BOO. iS ee eh ek 
_l 2000 


Ny; 320} 7 3 


240; 


160k + | 


” 


“Ny #100 
.Q +3000» 


* {d) 
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4.4.2 (Cont'd) 


asinapuremembrane, This effect is evident even for high pressure loading, (Figures 4, 4, 2~3d 

and 4d). Inaddition, the magnitudes of these tensile membrane stress resultants increase with in- 
creasing temperature difference and normal pressure while they decrease with increasing 
average temperature (quantities proportional to N -: This is to be expected, since increasing 


5° pressure and temperature difference each cause additional middle plane stretching while 
Nin tends to neutralize this effect. The membrane tension Nw decreases locally as the outer 


boundary is approached, where the gradients are most spnsenect: This negative increment 
AN.» may be responsible for the abrupt reduction in tensile hoop stresses (N,) in the vicinity 


of ee edge (Figure 4. 4, 2-4), For the larger values of average temperature, the high tensile stress 
resultants N in the plate interior reduce sharply in a boundary layer. Low temperature 


differences permit transition to compression near the boundary (Figures 4.4. 2~4b and d). 


(c) Bending Moments (Figures 4. 4,2~5 and ~6) 


The bending moments in the radial and tangential directions are constant 
and essentially equal over the major central area of the plate. However, radial moments 
decrease radically near the boundary, satisfying the zero moment boundary condition. Due 
to the predominantly flat profile of the deflected plate in the interior for large temperature 
differences (Figures 4, 4, 2-2a andb), it maybe conjectured that the constant and equal interior 
moments are the same as would occur in a clamped plate subjected to the same temperature 
gradient (since such a plate will remain flat), To show that this is the case, we proceed as 
follows: 


The bending moments in a fully clamped plate subjected to a thermal 
gradient through the thickness are given by (Reference 4-8) 


M.. 
MMe a 
while 
wz 6 


or, in nondimenstonal form, using the notation of (14) of Paragraph 4. 3,3 and (1) 


This result is readily verified by Figures 4.4.2-5 and -6. 
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SECTION 5 - AXISYMMETRIC STRESSES AND DEFLECTIONS IN SHELLS 


DUE TO THERMAL AND MECHANiCAL LOADS 


5.1 SUMMARY 


This section considers the axisymmetric stresses and defiections in heated snells 
under load. Basic equations, applicavle to general shells of revolution, are derived in 
Sub-section 5.3, These are discussed and applied specifically to conical shells (Sub-section 
5,4) and cylindrical shells (Sub-section 5,5) where solutions to the differential equations 
are given, 


0, 2, INTRODUCTION 


The purpose of this section isto present the theory fer the linear elastic analysis of 
shel's of revolution subjected to axisymmetric mechanical loads and temperatures, The 
formulation used is the linearized version of the development given by E, Reissner (Refer~ 
ence 5-1), extended to include the effects of temperature, 


The present work is exact withtu the framework of linear shell theory and removes 
the following restrictions of the more elementary presentation in Section 7, Volume I (Ref- 
erence 5-2): 


(1) The material in Volume Twas limited to the cases of truncated conical and 
cylindrical shells. 


(2) The development in that volume was approximate for the conical shelts and, 


(3) Thetemperature was constant through the thickness, i.e,, only meridional 
variation was permitted. 


Sinee the main objective of this section is to develop the seneral equakious governing 
the analysis of axisymmetric shells, nuinimat emphasis is placed upon obtaining specific 
numerieal results. However, problems invelving interaction between shells and bulkheads 
are discussed and, far the purposes of illustration, the general solution fer conical and 
cylindrical shells are derived, Numerical results are viven fer an unrestrained eylindri- 
cal shell subjected (o a prescribed temperature variation. 


The following symbols are used throughout this section: 


h Shell thickness 
a 
‘ ses Rh 


* 


. t=n 


ig Distance from a general point en the mertdian to the axis ef revelution 


8 Meridtonal coardinats 


Dae (Cont'd) 


u Displacement in the r direction 
Ww Axiai displacement 
z Axial coordinate 
3 
D Flexural rigidity, — th 
12(1-v*) 


E Young's modulus 

G(s), H{s) Load - pOnIpeTavUre functions defined by Eqs, (ib) of Paragraph 5, 5.2 and (1b) of 
Paragraph 5,4,2, respectively 

H Horizontal ferce per unit of circumferential length 

wR 


Curvature 
L Length 
M Moment per unit of length 
N Force per unit of length 
Maps Nop \. EATUA?, i EATdt, respectively 
p Surface tr nein (force per unit ef area) 
Pie PY Horizontal and vertical components c. P, respectively 
Q Shear foree per unit of circumferential length perpendicular to shell mid-plane 
R Radius of ¢ linder 
T Temperature cise above unstressed undeflected datum 
Vv Vertical forre per unit of circumferential length 
i ey » 

eX ig Gun Wale Og 
B Rotation of meridional element 
€ Strata 
t Thickness coardinate 
8 Azimuth or hoop angle 
d Coefficient of linear thermal expansion 
Vv Poisson's ratio 
£ Meriifonal coordinate 
A { san’ 
og stress 
QO Angle between tangent to meridian and the horizontal 

1 

2 
e eps 

SUBSCRIPTS 

o, 1 Refer te shell edges 
e, p Complementary and particular solution, respectively 
t, f Initial and final respectively 
8 In the meridional direction 
8 Tn the azimuth or hoop direction 
é In the meridional direction 
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5.3 BASIC EQUATIONS 


5.3.1 Configuration of a Typical Shell Element 


Since this study is concerned with axisymmetric problems for a shell of revolution it 
is sufficient to consider a typical meridional element before and after deformation due to 
loads and temperature (Figure 5,3.1-1). Positive directions for loads and displacements 
are as shown in the figure, The equation of the meridional curve is expressed parametricaily 
by 
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where £ is the independent variable which defines the position of a general point on the 
meridional mid-surface curve. However, tne quantity § may not have the units of length (for 
example, £ may denote an angle). Jn order to define length along the curve in a general man~ 
ner, we introduce the quantity @ = a(&), such that a differential arc length ds is given by 


ds = ade. (1b) 
But, 
pen = eee a 
eso ds ad& a 
(1e) 
ain a dz = dz = an 
PO ds ade 
Therefore, 
2 Brae 
(ry tz) oa, (1d) 


which defines a, 


We consider that the tangent to the meridionat curve befere deformation makes an 
angleo with a radial Une normal to the axis of revolution: the angle after deformation is 
designated by (07 8). Other quantities in the figure are self-explanatory or indicated in 
the nomenclature, and moment vectors are obtained usity, the right-hand rule, 


5.3.2 Equilibrium Fquattons 


In the following derivation the forces per unit of circumferential length are designated 
by the components (N .,Q) oF the statically equivalent components (V,H). ‘The relationships 


between these two sets of forces are given by 
Q=Vicos @- Hsing 
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where, inthe linear cheory, the initial configuration is emploved in writing the forec 
equilibrium equations (terms in & neglected). 


A free body diagram showing equilibrium of forces (Figure 5,3. 2-1) vields the fol- 
lowing equations: 
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FIGURE 5. 3.1-1 CONFIGURATION OF A TYPICAL MERIDIONAL ELEMENTS 
BEFORE AND AFTER DEFORMATION 
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FIGURE 5.3. 2-1 EQUILIBRIUM OF FORCES AND MOMENTS 


* This represents the vector sum of the hoop moments shown in Figure 5.3.1-1. Since the 
vectors representing the hoop moments form an angle cos od8 neglecting terma in 8), 
thetr resultant is of magnitude Mycos pdBadé with the direction as Indicated in this figure. 
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eVyyis rere so (2a) 
(vH)! ~ aN,t raP,, =0. (2b) 


From the free body diagram of Figure 5.3. 2-1b, moment equilibrium results in 


(rM,) - a cos” Mg ~arQ=0. (38) 


5.3.3 Curvature Changes and Mid-Surface Strains in Terms of Displacements 


Consider a meridional element of the shell mid~surface before deformation (Figure 
5,3, 1-1). The length of this element is given by 


=~ d 
dL, de 88 rag ‘ 
t cos o cos ~ 


After deformation, the new length is 


_(rf tut) d 
dL = cos (o + 5 


where 8 is the rotation of the element. Thus, the meridional strain of the mid-surface is 
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If we neglect terms higher than order 8, the linear meridional midsurface atrain becomes* 


a! sino 
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* Tt should be noted that the derivations In this section involve division by quantities r and r', 
Subsequent derivations also tnvolve division by z'. Therefore, at points where these quantl- 
tles become zero, the resulting derivations and equations become valid only when appropriate 
limiting processes and regularity conditions are imposed (see Paragraphs 5, 4 and §. 5), 
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The hoop strain of the mid-surface ¢g is a measure of the increase in circumferential 


length per unit of length and is therefore given by 


eh 
&g =r" (2) 


In order to obtain bending moments in terms of deformations, it is necessary to de~ 
termine the changes in principal curvatures of the shell mid-surface. From the definition of 
curvature, the change in meridional curvature due to deformaticn is 
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The hoop principal radius of curvature is defined by the length of the normal line tc 


the surface which ig bounded by the surface and the axis of revolution. Thus, the change in 
hoop curvature due to deform tion is 


K = sin @ + B) - sing 
8 r 


Again, neglecting terms of order higher than 6 yields 
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5,3.4 Stress-Strain and Moment-Curvatuwre Relations 


The stress-strain relations for the shell mid-surface including temperature are 


given by 
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where A is the coefficient of linear thermal expansion. Assuming that the strains are linear 
through the thickness of the shell, the bending strains are odd functions of Gand therefore in- 
tegration of the stress-strain reiations through the thickness ylelds 
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For purposes of obtaining expressions for the moments, it is further assumed that 
normals to the mid~surface before deformation remain normal after deformution, and tha 
mid~surface stretching due to bending is negligible. Then, employing the stress-strain 
law, integration of the first moments of the stresses through the thickness results in 


— =— a _T 
M, i o Gat D i Ky + 0K | es? 
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a3 
D = Eh _ 
W241 - v7) 
M. = | EXTYAC . 
h 
Substitutfon of Eqs. (1) - (3) of Paragraph 5.3.3 into Eqs. (2) and (3) gives the force re- 
sultants and moments in terms of slopes and displacements: 
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r , 
where Cos ~ = <= has bean employed. 


A compatibility condition in terras of force resultants and slope may now be obtained. 
Elimination of the displacoment u between Eqs, (1) and (2) of Paragraph 3,3 results in 
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Tha above nay be combined with Eqs. (2) to yield 
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5.3.5 Formulation of the Boundary Value Problem 


(4) Differential Equation 


In the following development, it will be assumed that V is known at one of the 
shell edges (V = Me at E = gO) so that from static equilibrium Eq, (2a) of Paragraph 5.3.2 


Tye Ad r@ PY dé. (1) 
eo 


This implies that ¥ is a mown quantity throughout the shell. 


Then the quantities @, Ne No. H, Me ze may be eliminated from among the 


seven equations (1), (23), (2b), (3) of Paragraph 5.3.2, and (4b) and (5b) of Paragraph 5, 3. 4, 
to vield a differential equation for the unknown slope 8. After detailed calculation and con~ 
siderable simplification, this equation may be written in the form: 
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where 


vinsage [Pe8(E) ee) -E) 


The Equation (2) is the general differential equation governing the linear analysis for stresses, 
and deformations in arbitrary shells of revolution due to axisymme' ‘ic loads and temperature 
(see, however, remarks made in footnote, page 5.7), 

(2) Force Resultant, Moment, and Displacement Boundary Conditions in Terms of 
the Slope 8. 

Since the differential equation (2) is in terms of the dependent variable 8, it is 

advantageous to express the force resultants, moments and displacements, as well as the 
boundary conditions, as functions of this variable. The derived quantities are given by 
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Typical boundary conditions can be expressed in terms of the above derived 
quantities as shown by the following examples: 


(a) Clamped Edge at & = os 
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In order to apply the differential equation (2) and boundary conditions GA toa 
specific shell geometry, the quantities r, @, z must be determined from the equation of the 
shell meridian, Specific examples for the cases of conical and eviindrical shells are pre- 
sented in the following paragraphs, In general, the maj t ,coblem encountered in solving 
Eq. (2) for a given shell of revolution involves obtaining the complementary solution for 2 
fourth order differential equation wits variable coelficients. Usually ciosed form solutions 
are not possible. However, as may be observed from ati inspection of the right hang side of 
Eq. (2), the introduction of temperature does not complicate matters, since it enters ina 
form andlogous to terms resulting from mechanical loads. 


Vhe boundary con fitions (Eq. (4d)) are of impertinee in considering the inter- 
action of the sheit with both transverse circular be'kheads wna other shells ef revolution. 
Once the linear response to edge load. and moments is determined from a solution of Eq. (2) 
subject to boundary conditions of the type (5d), intluence coefficients are known and internal 
loads at the Junctions of shells and bulkLeads may be evaluated by traposipg compatibility and 
equilibrium conditions on slopes and deflections. A description af the procedure to be fol- 
lowed is given in detail in Vohime I, Section &. 


o. CONICAL SHELUS 


o,4. 1 Basic Rqu.tions 


4 


A meridional section ef the cone and the coordinate syvsccm is shown in Figure S.-h 1-4, 


FIGURE o.d.d-l CONE GEOMETRY 


Referring fo Eqs, (1) of Parayraph 5.5.1, we choose the independent variable as 
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from which it follows that 
r =g§coso 
Z =s Sino 


r'=cos © (Lb) 
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where o> constant kas been employed and primed qiantiies indicate differentiation with re- 
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Correspondingly, ‘he expressions for the force resultants and mome {8 «re obtained 
from Eqs. (3) of Paragraph 5.3.5. «are given by 
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Then the homogeneous form of Eq. (1) becomes 
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where Jo and Yo are Bessel functions of the first and second kind of order 2. After much 


manipulation Eq, (4) may be written in terms of Kelvin functions as 
= meres 2 , ; i i+ 24° 
BA (rer j bord }* B eee "7 ber v) 


+ C ker ¥ - nai ¥)+ P (aut¥+ § ker (5) 


4 
where ° ~ hh 7 


The particular solution of Eq, (1) van be easily obtained if the function H(s) is ex~ 
pressed as a polynomial of the form 
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Then, selecting the solution on the form of series involving integral powers of 8, there 
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where the a, are given by 


P= A (p+l) (p) (m+) (m) 2 
_ 2 (p-1)! 
i a a p= (p-m+2) feat | (8) 


The total solution to the cone problem is 8 = B. + B,, and the constants A,B,C, F must 


be determined from the four boundary conditions (two at each edge). For the special case of 
a full cone (8, = 0), finiteness of stresses at the apex requires that La = 0, Further, regu- 


larity conditions on stress resultants and shears, require that "C' ="F" = 0. Specific 
solutions for the full cone then resolve themselves into the determination of the two con- 
stants ''A" and ''B" from the boundary conditions at the base. 


The Kelvin functions and their first derivatives appearing in Eq. (5) are extensively 
aud accurately tabulated in Reference 5-4 for a wide range of the argument. The form of 
solution developed above ig readily adaptabie to numerical computation using digital com- 
puters. Parametric stlidies to determine stresses and deflections can be made. 


5.5 CYLINDRICAL SHELLS 
5.9.1 Basic Equations 


The cylinder and the coordinate system is shown in Figure 5, 5. 1-1, 


FIGURE 5.5,1-1 CYLINDRICAL GROMETRY 


We choose 8 again ag the imependent variable , however, since r' = 0 for the 
cylinder (see Figure 5. 3,1~1), then as indicated previously Eq, (2) of Paragraph 5.3.5 {s 
not applicable in the form given. The equations for the cylindrical shell are more readily 
derivable from the conical shell equations (which were developed from Eq, (2) of Para- 
graph 5.3.5), Proceeding in this manner, {f the equation (3) of Paragraph 5.4.1 ia multi- 
plied by 2, 
taken, there results, for the caso of normal pressures, 


the quantity + replaced by aoe. and the appropriate I{mit as orn, ro R, {s 


t 
M,,.' N,, 
pal” + EM p- — Pa f (1) 
R 


5.5.1 (Cont'd) 


Similarly, the limiting process applied to Eqs. (4) of Paragraph 5.4.1 for the force 
and moment resultants yields 
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An alteruate form of the equilibrium equation in terms of the racial displacement com- 
yonent "uu can be obtained by integrating Fq. (1) once, noting, that 8 = -ut, This results in 
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where the constant of integration is giver by the last term on the right. Substituting B= -u' into 
Eqs. (2) the force and moment resultants are 
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Using the formulation ot Eq (1), for all sets of boundary conditions involving B, the 
solution may be integrated once to obtain u. The constant of integration is determined from 
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the prescribed value of Vo When using the formulation of Eq, (3), however, a appears 


explicitly in the differential equation, 
5.5.2 Solution of the Differential Equatior 


Equation (3) may be written as 
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A particular solution can easily be obtained if G(s) is expressed as a polynomial with 
terms of the form 
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Then a particular solution corresponding to G{s} = o£ 


this solution in the form of a polynomial in (8/2) of degree K and substituting in Eq. (1). 
There results 
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The constants Ai A,» A,, A, must be determined from the specified boundary conditions at 
he i 
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the ends of the cylinder, This requires the solution of four simultaneous, linear algebraic 
equations, In order to eliminute the necessity of solving four equations an approximate method 
of solution in which the less tedious procedure of solving two pairs of simultaneous equations, 
pach for two unknown constants, is used extensively, Tne basis of this method is the assump- 
tion that edwe shears and moments (which are self equilibrating, applied at one end of the shell 
negligibly affect the stresses and deflections at the opposite end, This assumption is valid if 


the length of the shell 'L'' satisfies the condition I> 3/, or equivalently, for v- 30, when 
” 


ee aoe Beare page cy es ‘ te yee 
¥ * < ce) . The inequality is satisfied in all but very short/thick cylindrical shells. 
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The method, which appears in Referenee 5-2, may be described as follows: 


(i) Determine the detlections and stresses corresponding to the particular solution 
M 
. 
on te ) (u, \ where the{ uy) are given by Eq, (4) In general this solution will not satisfy 
k=O : a 


the boundary conditions. 


2) Edge maments and loads are then applied to the semi-infinite cylinder such that 
when their effects are superposed on the solutions of (1) the boundary conditiona at each edge 
are satistied, Since ibere are only two conditions to be met at each edge and there is no assumed 
interaction of edge effects then, for each edze, the solution of two simultaneous equations in two 
unknowns (one shear and one bending moment) is required, 


5.56.3 NUMERICAL EXAMPLE 


A free cylindrical shell is subjected to an elevated temperature as shown in Figure 
5.0,3-1, ia, the inside surface of the cylinder is ata uniform clevated temperature T, and 


the outside temperature varies linearly from T, at s*) to a temperature T. (itA) at s-L. The 
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variation through the thickness of the sheil is assumed to be linear. It is required to de- 
termine the stresses and deflections. 


FIGURE 5.5.3-1 ILLUSTRATIVE EXAMPLE 


Proeedure: 


The temperature distribution is given by 


T(s, 0) = 4, E Ae (- 3] 


and therefore 
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Then from Eqs. (2), (43, and (5) of Paragraph 5. 5.2 the solution for the displacement is 
given by 
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The boundary conditions for free edges are 


M os Q SO at os. oO, hay 
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whieh from Faqs, 4 of Puragraph 5.5, 1 can be written as 
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The constants Ay - am in Eq. (1) are determined from the boundary ev cditions. 


Foree resultanc. and moments are then found by substituting into Eqs. (2) of Paragraph 
5.5.1, Sinee the caleulation details are straightforward, only the results are shown, 
Figure 5,5, 3-2 vives nondimensional deflections, force and moment resultants for both 
acylinder with Lf ° 10 and a ionger cylinder corresponding te L/¢ = 50. The graphs 
show that as the cylinder becomes longer, the peak deflections and stresses approach the 
end of the surface subjected to the higher thermal gradient. In general, these peak values 
tend ta increase in magnitude sor the longer eylinders, resulting in sharp gradients in the 


vicinity ef the edge. 
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